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»' SIMULATION AND OPTIMIZATION OP A WATER-GAS 

SHIPT REACTOR” 

"by 

K.K. PAHUJA 

Department of Chemical Engineering 
Indian Institute of Technology, Kanpur 

ABSTR/iCT 

A simulation model for the water-gas shift reactor 
is developed using the kinetic data for the PCI H. T. shift 

f; 

catalyst to estimate the values of the reaction rate constant 
The model takes into accoxint the effect of diffusion within I 
the catalyst pores. The model appears capable of supplying 
resxilts in reasonable agreement with the plant measurements. | 
The model is used next to prepare an optimization | 

scheme for optimizing an existing industrial reactor with | 

1 

maximum CO-conversion as the criterion. The computations 
have been made on IBM 7044 computer. An improvement in the 
plant conversion from 87. 97?^ to 91.14^ is predicted at the 
given feed composition and the reactor load. The effect of 
the various parameters like steam/gas ratio, dry gas feed ' 
rate and catalyst age on the optimum conversion is also studie 
Over the range studied, the optimum CO-conversion increases 
with increasing steam/gas ratio, decreases linearly with 
increasing dry gas feed rate and is hardly affected by cataly 
ageing beyond a period of 270 days. 





CHAPTER 1 


INTRODUCTION 

Th© recenl; "trGnd in the study of process systems is 
througH simulation. Simulation can be broadly defined as 
representing a physical process or phenomenon by means of a 
mathematical model. 

Off-line computer simulation of an actual plant can 
provide a priori knowledge of the behaviour of the plant, 
provided that the mathematical model chosen for simulation 
describes the plant with reasonable accuracy. The simulation 
model achieved by writing down the process equations based on 
the physical and chemical operations involved in the plant is 
programmed on a relatively large off-line digital computer 
and experiments are performed on the simulated plant to observe 
the effect of changing independent variables on the plant 
performance. This type of computer simulation is faster and 
cheapter than experimentation on a running plant. In addition, 
the range of independent variables for simulation can be larger 
than that observed normally in a plant in a day-to-day operation 

In the present work, a simulation study on an industrial 
shift reactor has been carried out with the purpose of using 
the simulation model in conjunction with an optimization scheme 
for determining the key variables that will yield maximum 
conversion. 
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Description of the Reactor ; 

The industrial reactor concerned in the present work 
is the multi-hed CO converter of the fertiliaer Corporation 
of India, Tromhay. The reactor is shown schematically in 
Fig, 1,1, In the preparation of ammonia synthesis gas, after 
partial oxidation of Naphtha, the product gas contains hydrogen 
and carbon monoxide as the major components. In the converter, 
this carbon monoxide is converted with steam to produce 
additional hydrogen via the shift reaction 

CO + H20(g) - COg + 

The reaction is carried out in the temperature range 
of 330--500°C in the presence of Chromia promoted iron-oxide 
catalyst (high temperature shift catalyst) developed by 
Fertilizer Corporation of India - Planning and Development 
Division, 

Before entering the converter, the gas is saturated 
with water in a heater- saturator, mixed with steam and is 
further heated in a preheater. The converter consists of 
three catalyst beds, which are essentially adiabatic. Usually 
the first stage is operated at a higher temperature to get 
fast rate of reaction, and the last stage is operated at a 
lower temperature to increase the equilibrium conversion of GO. 
For removing the exothermic heat of reaction, the gas between 
the beds is cooled in interstage coolers. 






4 


CHAPTEH 2 
THEORY 

I, Steps in Gas-Solid Catalytic Reactions : 

The ma 3 or role played hy the catalyst In a reaction 
is that it increases the rate of reaction because it makes 
possible an alternative mechanism, each step of which has a 
lower free energy of activation than that for the uncatalyzed 
process, without affecting the state of equilibrium between the 
products and reactants. 

In a gaseous reaction within a packed-bed catalytic 
reactor, the various physical and chemical processes have to 
be taken into account while considering the overall catalytic 
process. The different steps involved are: 

1 , Mass- transfer of reactants from the fluid phase 

to the external surface of the catalyst particles. 

2, Diffusion of reactants into the pores of the 
catalyst particles, 

3, Adsorption of reactants on the catalyst surface, 

4, Reaction on the surface, 

5, Desorption of products from the catalyst particles, 

6, Diffusion of products out from the pores of the 
catalyst particles, 

7, Mass-transfer of products from the external surface 
into the fluid phase. 
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_ The steps 1, 2, 6, and. 7 are entirely of a physical 
nature, whereas the steps 3, 4 and 5 are chemical phenomenon^ 

External Diffusion ; 

The rate of diffusion, r^, can he expressed by the 

relation 

rj, = kg a (Pg - p.) ( 2 . 1 ) 

Where 

Pg = partial pressure of diffusing gas in bulk stream 

% ~ partial pressure of diffusing gas at solid surface 

k = diffusion rate constant or coefficient 
s • 

a ~ external area of catalyst particle per unit mass. 

The diffusion-rate constant depends upon the nature i 
of the diffusing component, turbulence conditions in the vicinity 
of the catalyst particle and properties of the entire gas mixture 
The effect of bulk diffusion is the most important when rate of : 
reaction is high and linear gas velocity is low. In fixed- I 
bed reactors, where the fluid velocity is high with respect 
to the solid, molecular diffusion will play a small role in i 

determining the rate of mass transfer, 

III, Internal Diffusion in Porous Catalysts ; I 

Highly active catalyst pellets obtained by compressing 
finely divided powder, are porous and usually possess a large 
amount of internal surface. It is necessary for the reagents 
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to succeed in reacting this internal surface, and this will 
involve diffusion within the pores. 

The measured rate of reaction on a porous catalyst 
includes the resistance of internal diffusion. The factors 
which increase the importance of diffusion are an active 
catalyst surface, a high temperature, a large particle size, 
and small pore diameter. When the diffusion resistance is 
significant, the activation energy, form of the rate equation 
(order of reaction), and other kinetic factors become depen- 
dent upon pore size and particle size. The internal diffusion 
results in an apparent decrease of the energy of activation 
with increasing temperature, 

» 

Two kinds of diffusion in pores are possible. If the 
pores are large and the gas density is high, the process is 
that of bulk, or ordinary diffusion. The effective diffusion 
coefficient, will be given by 

^eff ~ /T (2.2) 

¥/here B is the ordinary diffusion coefficient, S is the 
voidage and T is the tortuosity factor which has to be intro- 
duced due to the pores being neither strai^t nor of uniform 
cross-section. The diffusivity is inversely proi)ortional to the 
pressure, 

If the pores are quite small, and/or the gas density 
is low, the molecules collide with the pore wall much more 
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frcquBn tly ^ tlian with each other* This is known as ^ Knud sen 
diffusion*# The Knudsen diffusion coefficient for a straight 
round pore is given by 

Djj. = 9700 rg /t/M CmVsec (2.3) 

1 

Where 

Tq = equivalent radius of pore, era. 

T = Temperature, °K 

M = Molecular weight of the component. 

:i^,l 0^ this diffusivity, 

IV. Catalyst Effectiveness ? 

Since diffusion requires a concentration gradient, 
the concentration of reactants will be less in the centre of a 
catalyst particle than at its outer surface. Hence the 
reaction rate will in general be less at the centre than at . 
the outside surface. 

A quantitative evaluation of the influence of pore 
diffusion on the rate of reaction can be performed by introducing 
the term 'effectiveness factor', E. This is defined as the 
ratio of the observed rate to that which would be obtained 
if the whole of the internal surface of the pellet were available 
to the reagents at the sane concentration as they have at the 
external surface. 
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The mathematical analysis for determining the effective- 
ness of a porous catalyst in case of simultaneous diffusion 
and reaction becomes simplified by considering a spherical 
catalyst pellet and assuming isothermal conditions. It is 
further assumed that the complicated diffusion phenomenon 
within the porous structure can be represented by a single 
overall effective diffusion coefficient, D 

0XX • 

hor a reaction in which there is no change in number 
of moles on reaction, a steady state mass balance on a differen- 
tial spherical shell of radius IL and thickness dl? gives 


d^C ^ _ r 

dE^ ®p ^eff 


(2.4) 


V/here r = reaction rate per unit volume of catalyst pellet. 

C = concentration of reagent. 

The above eq,uo.tion is to be solved for the boundary conditions; 
1) 0 = O3 at Hp = Eg 

(i.e. the outside surface of the pellet). 


iU 11^ = 0 at Ep = 0 ; 

(i.o, no concentration gradient will exist at the ! 

centre of the sphere). i 

Solution of the differential eqn. (2,4) subject to the 
given boundary conditions describes the concentration profile 
in the catalyst pellet. For reactions of other than first i 

order, analytical solution of this equation becomes more difficult.! 
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V, The Qptimigation of Exothermic Equilibriutn Reactions : 

The optimization of the reactor output - i,e, maxi- 
mization of the amount of reaction product per unit of time 

and reactor volume may play an important part in case of 

(1 \ ’ 

reversible exothermic reactions \ Por such reactions the 
effect of a temperature increase is to raise the speed of the 
forward reaction but also to lower the value of the equilibrium 
constant and thereby to reduce the maximum attainable conversion. 
Thus the kinetic and thermodynamic factors are, in a certain 
sense, in opposition, A compromise must then be found between 
the requirements of a high degree of conversion (low temperature) 
and those of a high average reaction rate. Near the inlet to 
the tubular reactor, where the reacting gas is still far from 
equilibrium, it is advantageous to use a fairly high tempera- 
ture in order to increase the forward reaction rate and thereby 
to raise the reactor output. Near the outlet to the reactor, 
however, the temperature should be brought to a much lower^ 
level so that a better equilibrium state may be approached. 

Por a reversible exothermic reaction, the rate will 
pass through a maximum as a function of temperature when the 
composition of the reaction mixture is constant_ (Pig. 2.1) . 

At a low temperature, the reaction 'rate, r, will be 
small, whereas a high temperature, T*, will exist at which the 
mixture is in equilibrium and the rate is zero. The temperature 
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at wliich the rate is a maximum will be somewhat lower than T* 
and is determined by the condition 

\ ■ 
tf = 0 _ (2.5) 

The optimum temperature profile in a tubular reactor 
can be found on the assumption that the criterion (2,5) must 
be valid for any cross section of the reactor. Then the 
reaction rate has its maximum value as a function of tempera- 
ture at any degree of conversion. However in practice, the 
establishment of the optimum temperature profile to maximize 
the reactor output, involves technical difficulties and 
expensive provisions because of the very special schedule of 
heat removal which would be required, For gas reactions in 
a catalytic reactor, it is much simpler from a technical point 
of viev/, to let the reaction proceed adiabatically . A reasona- 
bly close approximation to the maximum rate path can then be 
achieved by carrying out the reaction in a multistage tubular 
reactor with interstage cooling of the reaction mixture by 
means of heat exchangers (Eef. Fig. 2,2), 

Each catalyst section operates adiabatically and 
consequently the temperature rises from inlet to outlet which, 
of course, is contrary to the theoretical ideal-. However, 
by removing heat from the gas after it leaves each section, 
the temperature Is brought to a much lower level at the entrance 
to the next section in the series. Thus an approximation to 




Fig. 2-2 - Multi-stage adiabatic tubular reactor with cooling 
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the optimum sequence is obtained by stepwise change. HOM 

( 2 ) “ 
and KUCHLER' provided the solution to the problem of opti- 
mizing such a system, i,e, that. of determining the highest 
possible final conversion for a given number of sections and 
for a given entrance temperature of the first section. (This 
is equivalent to determining the temperature sequence which 
will result in a minimum quantity of catalyst, for the attain- 

m * 

ment of some chosen hourly outpt of product^^^). The inlet 
temperature of the first stage, T^^q, is an independent para- 
meter which has to be varied in order to find the minimum 
amount of catalyst needed for a certain final degree of 
conversion, or the maximum degree of conversion obtainable with 
a given amount of catalyst in a given number of sections, Eor 
an U-bed sequence (Pig. 2. 2) the conditions for this optimization 
as given by Kramer^ are (Ref, Appendix I) s 
i’)' Por the optimum temperature policy, the condition for 
each section (index n) is; 



paths. 

ii) Between two sections (n-1 ) and n; 
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^(n-1),L = ^n,0 n = 2,3 ,IJ 

(2.7) 

Ihis condition relates to t since for the inter- 

stage coolers it is clear that 

^(n-l),L "" ^n,0 (2*7a) 

and 

^n,0 -<1 ^(n-1),L‘ (2.7b) 

The criteria represented by equation (2,7) implies 
that the rate of reaction at the exit of a bed should be the 
same as that at the inlet to the next. This is possible because 
of the double valued nature of the rate function (Ref ,T'ig.2, 1 ) . 

It may occur that T^ exceeds the maximum permissible 
temperature, so that only a corresponding degree of conversion 
(-C'x.l j^) cjin be reached in the first section. In that case, 
eqn,(2,7) vail not apply any more for finding the proper 
entrance temperature to the second section. 
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LITERATURE SURVEY 

Chemical Reaction Rate of the Shift Reaction ; 

A number of workers have studied the behaviour of 
the water-gas shift reaction over high temperature iron- based 
shift catalyst. Various forms of empirical rate equations 
have been suggested by different authors. A few investigators 
have proposed and tested reaction models also. 

Laupichler^^^ used a simple bimolecular expression for 
determining the rate of the shift reaction in the homogeneous 
phase. In the process of an industrial-scale production of 
hydrogen, in which a large excess of cteam is added to the dry 
gas, the change in concentration of water vapor can be neglected 
during the course of reaction, Eased on this fact, Laupichler 
justified the use of a pseudo first order rate equation in 
carbon monoxide in place of second order rate equation. 

It was considered that the effect of the catalyst was 
to increase the rate of conversion without changing the form 
of the rate expression. The effect of external diffusion was 
also taken into consideration, 

Padovani and Lotteri^^^ also assumed a second order 
expression for the forward rate of reaction, r^, i.e. 

^f = ^ Pqo 


( 3 . 1 ) 
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Where 

k = reaction rate constant 
Pj[_ = partial pressure of the component i 

Kulkova and Temkin^'^^ found that the following rate 
expression adequately covered the experimental data. 

r = k Poo ( I'OOj'Sh/I’HjoP''’ 

(3.2) 

with, h = 0.44 to 0,50 

According to this equation, H 2 has a retarding effect 
on the rate of forward reaction, whereas CO 2 has no such 
influence. 

Atwood^ and Mars^^^ applied a first order kinetic 
expression of the following type; 

r = k (Pqq - p5q ) (3.5) 

Where 

p*Q = partial pressure of CO at equilibrium. 

Popov^"**^^ used an unpromoted iron, oxide catalyst and 
found that the rate was first order with respect to CO and 
zero order with respect to H^O. At temperatures above 400 ^ 0 , 
the reaction order with respect to H 2 O changed to 0.5. 

( 11 ^ 

Barkley and coworkers^ ^ studied the kinetics of the 
reverse shift reaction using the - CuO shift catalyst 

and obtained the following rate equation: 
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1 + Kqq Pqq + Pqq^ 


It is not very clear if this equation will hold for 
the water-gas shift reaction also* 

( a ?'] 

Kodama^ ^ studied reaction rates over a coprecipitated 
^^2*^3 ” ^^2*^3 arrived at an equation of the 

type: 


Poo ^ H20 " ^002 %2 


(3.5) 


In this study, the contact times were far shorter 
than those used in industrial reactors. 

Later Kodama^^^^ using an industrial catalyst of some- 
what different composition arrived at the following rate 
expression: I 

PCO ^2^ ^ ^^2 I 

^ ^ 1 + K(30 Pqo %2° ^'^2° ^^2 ^^2 ^002 P 0 O 2 i 

(3.6) I 

■ 

The latter equation was found to he valid over larger j 

concentration ranges than the former. I 

Bortolini'' studied the shift reaction at atmospheric 

pressure and at relatively high temperatures using a commercial! 

Pe„0^ - Gr«0, catalyst. It was found that at temperatures of, I 
2 3 2 3 • ! 
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about 500 °C and higher, the reaction rate is controlled by 
the external diffusion. 

Stelling and Krusenstierna^/^^^ investigated, the 
reaction rate in the presence of an iron-oxide catalyst. It 
was found that even very small amounts of OO2 have a strong 
inhibiting effect on the forward rate whereas H2 does not 
take part in the reaction. Their rate equation for the 
forward reaction is: 


r^ = k 


^00 % 2 ° 


PCO2 ^ 


H2O PH2O 


(3.7) 


X16) 


Kirillov'’ ^ used an industrial ^620^ catalyst and 
developed the following rate equation; 


r 


Pen V 

^ Ph 0 ( ^ Ph ^ V — ^ ^ 

^2^ PcOo ^^2 Pco 


2 n'2 


(3.8) 

Atroshchenko and Bibr^'*'^^ studied the reaction rate 
over a commercial catalyst containing MgO and Pe20^ as the 
major constituents, in the pressure range of 6 to 21 atm. 
The data was correlated by the equation; 

Vo 


k 


Pco 


• * 

Pco 




^PH20^ 


(3.9) 


Hulburt and Srinivasan 


(18) 


used a commercial promoted 


iron oxide catalyst and arrived at the following equation for 
the forward rate at 350*^0: 
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0.2107 Phq 
% 0 

1+1,063 


( 3 . 10 ) 


According to this 
affected by and 


expression, reaction rate is not 
whereas it increases vath p^ . 


Bohlbro^'^^^ investigated the kinetics of the shift 
reaction at atmospheric pressure over an iron-oxide chtomium 
oxide catalyst, in the temperature range of 330-500®C. It 
was observed that the rate is independent of the hydrogen 
partial pressure, whereas an increase of carbon dioxide partial 
pressure serves to decrease the reaction rate. He proposed 
an empirical rate equation of the exponential type for the 
forward reaction. 


^f 


kCPco) 




The exponents were found to vary within certain ranges 
of temperature. The following ranges were obtained for the 
exponents; 

1 = 0.80 to 1.00 
m = 0. 20 to 0. 35 
n =-0.65 to -0.50 
q. = 0 

(' 20 ') 

Atroshchenko and coworkers^ ' also used the same 
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exponential type rate equation. The investigation was 
performed at 300°C and 460°G over an iron-chromium oxide 
catalyst, and different values for the exponents were obtained 

at the two temperatures. 

( 2 ^) 

Later Goodridge^ ' also investigated the same expo- 
nential rate expression using several iron base shift catalysts, 
in the temperature range of 350-430° 0. The exponential 
constants varied significantly as a function of temperature 
as well as with each of the catalyst. His study was, however, 
in agreement with Bohlbro's work for the same catalyst. 

found that Laupichler's second order 
expression was applicable for calculation of industrial 
converters using the iron oxide shift catalysts with 
as a promoter. The rate equation can be written as; 


r = k. 




K 



(3.12) 


Where, 

y^ = mole fraction of component i 

k^ = second order rate constant for the forward reaction 
K = Equilibrium constant for water gas shift 
reaction 

The rate constant was further modified to consider the 
effect of factors like size of the catalyst particles, period 
of run and operating pressure. 
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I’rom the investigations of different workers, it 
is evident that even with extensive kinetic studies for 
the shift reaction, there seems to be little agreement 
about the exact form of the rate equation. 



CHAPTER 4 


THE SIMLATION MODEL 

The simulation model is developed hy taking into 
consideration all the physical and chemical processes 
occurring in the reactor based on certain simplifying assum- 
ptions. 

The following assumptions are made in deriving the 

model. 

1. The bed is completely adiabatic. 

2. The composition and temperature in the 
radial direction are uniform. Because of 

the adiabatic nature of the bed, this assumption 
is justified. 

3. The velocity profile is flat at any cross 
section, A large ratio of bed diameter to 
the catalyst pellet diameter implies that the 
velocity profile ±b essentially uniform along 
the bed diameter, 

4. There is negligible axial diffusion in 
the bed. 

5. The reactor is in a steady-state. 

6. Pressure is constant over the length of 
the reactor. 

The assumptions 2, 3 and 4 constitute the plug-flow 



reactor model. 


I . Material Balance ; 

Assuming CO as the guide- species, a material balance 
over a cross section of the catalyst bed leads to 

^f,C0 • ^ ” (iV (4.1) 

Wie re , 

P = Total flow of gas mixture, Nm^/hr. 

j.y . = Mole fraction of component i in feed stream, 
i. , 1 

X = Degree of conversion of CO. 

r = Overall rate of reaction of CO , 

(m^ bed volume) (hr) 

5 

V = Volume of reactor bed, m . 


II, Energy Balance : 


The energy balance in case of adiabatic operation is 
obtained by equating the heat of reaction term to the sum 
of the sensible heat gain for the gases. Por a system with 
no change in number of moles on reaction, this yields 


i ^ 


c (T,P) 

Pi 


dT 



AHj^(T,P)J 


y 


f,CO 


• dx 


( 4 . 2 ) 


¥here , 


0 


P. 


= mole fraction of component i 

=: Heat of reaction, cal/gmmole of CO 

= molal heat capacity of component i, cal/(gmmole)- 
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*■ T = Absolute temperature, °K 

P = Total pressure, atm. 

The summation in the energy balance equation is over 
all the reactants and products present in the reaction mixture. 

The nonideal behaviour of the gases has been considered 
and, therefore, the values of heat capacities and heat 

of reaction 41 are dependent on temperature and pressure. 

The specific heats are expressed by Hougen and Watson^^^^ as 
polynomials in T; and their pressure dependence has been 
obtained by using data^^'^^ at different pressures (details 
are shown in Appendix II), Table 4.1 lists the heat capacity 
relationships for various gases; the pressure dependency for 
gases other than steam is negligible. 


1 . 

2 . 

3. 

4. 


5. 

6 . 


7. 


TABLE 4. 1 : RELATIOISHCPS POR HEAT CAPACITIES OP 
VARIOUS (LISES ^ ^ 

( T in °K ; C^ in g-cal/(gmmole )( ®K) ) 

Gas' ' Heat Capacity Equation* ' " ' ' 

H 2 X 10"^T +.0.4757 X 10"^ T^ 

CO (Cp )2 =6.350 + 1.811 x 10“^T - 0. 2675 x 10"^ T^ 

GO 2 (S^3 = -t'10.14 X 10“% f 3.415 X 10"® T^ 

(Cp)^ = 3.204 +18.41 X 10~% - 4.48 x 10“^ T^ 

% (Op)^ = 6,457 + 1..389 X 10“% - 0.069 x 10"^ T^ 

A (Op)g = 4.97 

HgO (^ 2^7 = 7, 136 + 2.640 x 10“% + 0.0459 x 10“^ T^ 

+ 0.63 X 10 "^^ p1.29yp5.56 

*Pressure range 1 to 40 atm. 
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The heat of reaction has also been expressed as a ! 

function of temperature and pressure (Details shown in 
Appendix II) and the resulting expression is 

=[Aa.T + ^ T^ + T^ + 13.82*10'’^ p1 . 29^-4-. 56^ 

- 2.96 - 0 . 933 * 10 '^ exp(-0.9354*10*’^P)-672.47 

(4.3) 

With 

Aa = -0.201 , ! 

Ab = 5.493 * 10~^ 

Ac = -2.7177 * 10“^ 

The above equation is valid in the normal operating 
range of shift reaction (Temperature range 600 to 800 ®K and 
pressure range 1 to 40 atm.). 

III. Kinetics of the Shift Reaction 

As appears from the literature survey (Chapter 3), 
a generally accepted kinetic expression for the shift reaction 
does not exist. It is quite justifiable, therefore, to make 
use of a kinetic equation vdiich is both simple and sufficiently ^ 
accurate at least for industrial purposes. Thus the rate 
equation chosen for developing the simulation model is the 
second order rate expression^^^ for which a pseudo first 
.order equation is a reasonable approximation when steam is 
present in considerable excess. The derivation is shown in 
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AppeMix III* Ihe pseud©-flrst orde^ rate expl^SiSion dan 
be written as 

^ ( ^ 00 " 5^00 ) 

with 


“ ^2 

Where 


(4.4a) 


y^Q-s mole fraction of CO 
y*QQ = equilibrium mole, fraction of CO 

“ second order and pseudo first order rate constants 


^HoO 


respectively, . 

= mean value for the mole fraction of H«0(g). 


(It will often be satisfactory to assume that a givbn bed 
operates at the mean of the inlet and outlet steam concentration 
However, for more accurate calculations, a stepwise procedure 
may be adopted where the mean concentration is estimated over 
small increments of reactor bed). 


C 25 ) 

It has been shown by Ruthven^ ^ that undd3? normal 
operating conditions, external mass transfer resistance is 
Tjaually negligible in an industrial shift converter,^ It 
appears' that under most conditions the major resistance, is 
due ■ t'o‘ pore diffusion of carbon monoxide. The assumpti 9 n of 
thermal equilibrium between the catalyst and gas with a uniform 
temperature throughout the catalyst pellet appears to be 
justified^ 
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The pore diffusion tends to alter the value of the 
effective rate of reaction and this modification can he taken 
into account by introducing a suitable effectiveness factor, 

E. For the direct evaluation of E in case of (pseudo) first 
order reaction, the correlations proposed by Thiele^ can 
be used. According to him, in case of a first order reaction 
taking place in an isothermal catalyst particle, the influence 
of a restricted internal diffusion on the reaction rate is 
governed by a dimensionless modulus m. Using the formulation 
of Hoogschagen'- the dimensionless modulus has the form 



in which 




catalyst particle diameter, cm, 
s intrinsic catalyst activity, sec ’ 

S3 effective diffusion coefficient in the porous 

p 

catalyst under reaction conditiotis,,.cm /sec. 


disregarding the change of with conversion, the value " ; 

of m is independent- of concentration for a first order reaction. 

From the solution of the diffusional equation (2,4), 
the effectiveness factor, E, of the catalyst pellet can be 
determined and this leads to the relation 


3 ( 1 . 1 ) 

m ^ tanh m m ^ 


(4.6) 
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The effective particle diameter, Dp, can he calculated from 
the formula 

Dp = 6 Yp/Ap (4.7) 

in which Ap is the gross exterior surface area and Vp is the 
volume. Por spheres this formula gives the actual diameter. 

The apparent reaction rate constant, ko^, based on 
unit catalyst bed volume is ordinarily different from kv^ 
and can be expressed as (29,26) 

ko^ = kv^ (1-e) B (4.8) 

Here e is the voidage of catalyst bed and its value can be 
taken as 0,4 for uniformly size 2 particles 

Por computing the apparent activity, ka^, defined 
with respect to space velocity measured under standard 
conditions, usxially N. T,P,, the relation used is^^^^ 

ka^ = ko^. 3600^ (^).? (4.9) 

The units of ko^ and ka^ in this equation being Sec”’^ and 
hr respectively. 

The question of the effective diffusion coefficient 
in ferric oxide-chromic oxide shift catalysts has been discussed 
by Hoogsehagen^^®^. It has been shown that for the catalyst 
tablets prepared by compression of a powder, the diffusion 
through the highly heterogeneous pore system predominantly has 
the character of normal gas diffusion. The value of effective 
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diffusion coefficient, was estimated to 0,.039 sq, 

cm/sec. for CO at 1 atm« and 350'’G„ It can be assumed to ■ 
a good approximation that Dg^j^ does not change vdth the 
degree of conversion over the range of gas compositions 
encountered in practical shift converters^ Assuming 
the same effective diffusivity data to be valid for the 
present catalyst, the effective diffusivity at other pressures 

/ on pg 

and temperatures may then be calculated from the expression' 


“eff = ( SB omysec (4.10) 

The value of from this equation can be used 

for calculating the Thiele modulus and therefore the effective- 
ness factor. 

The values of the kinetic parameters have been evaluated 
using the pilot plant kinetic-data reported by Chandra etal^^^^ 
for the FCI H.T. Shift ■ Catalyst, The pilot plant data have 
been obtained using 6 mm dia, * 3 mm. long catalyst tablets. 

Use of equation (4.7) gives = 0.45 cm. for which the pore 
diffusional resistance must be significant. 

When pore diffusional resistance is important, Eqn,(4. 6) 
simplifies to give 


Eqns, 


E Cii 1 (1 - 1 ) (4. 11) 

Valid for m >■ <^3, 

(4.8), (4.9) and (4. 11) can be combined to give 
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,p- 2D' 

^ ( J kv^. - — ^££)(i-e). 3600 (^)P (4.12) 

Eqns. (4.10) and (4.12) can then he used for evaluating 
the value of kv^ at a given pressure and temperature from 
the value of ka^ which is available from pilot plant data. 

Calculation of Intrinsic Activation Energy : 

Ihe temperature dependence of the intrinsic rate 
constant, kv^, can be expressed by a simple Arrhenius 
relationship, 

-E /RT 

kv^ = e (4.13) 

Where, 

E^ = intrinsic activation energy, gm-cal/gm-mole 
R = gas constant = 1.987 gm-cal/gm-mole °K 
T = Temperature, °K 

Z = Pre exponential factor of intrinsic rate 

' -1 
constant, sec 

The pilot plant data at two different catalyst life- 
periods and pressures have been used for determining the value 
of E„. The results calculated are shown in Tables 4.2 and 4.3. 

■a* 

By plotting log kv^ versus l/T, straight lines are 
obtained as shown in Big, 4.1 and 4.2, Erom the slope of the 
straight line, the value of obtained is 9.33 * 10^ and 

9,12 M- 10^ respectivelyo which are practically identical. 
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Fpoffl tJiese the average intrinsic activation energy hai Ijeen 
calculate^ to be 18.330 * 10^ gm ca?^/gm-mole, 

This value for the intrinsic activation energy 
corresponde fairly well with the latest values of 21.4*10^ 
to 27.3*10^ cal/iaole calculated by Ruthven^^^^ for various 
commercial shift catalysts^ 


TABLE 4.2; CALCULATION OF INTRINSIC CATALYST 

ACTlfxTlT'B? VARIOtfS feSASlfloM IPJEUBRATURES 


No.. 

Temperature 

OC 

Avg, value 
of experi- 
mentally 
determined 
rate cons- 
tant ka^, 

hr-''. 

TW 

Calculaied 
values of in- 
trinsic rate 
ccmstant^ 
kv^# se§-‘^ 

1. 

480 

0.3235 * 10^ 

1^328 * 10“^ 

8,36 

2^ 

460 

0.25965* 10^ 

1.364 * 10"^ 

5.48 

3. 

440 

0,2215 * 10^ 

1,403 * 10“^ 

4.05 

4. 

435 

0.2340 * 10^1 

1.412 * 10*^ 

4,435 

5. 

420 

0.18517* lO^i 

1.443 * 10“^ 

2,88 

6. 

415 

0.1670 * 10^ 

1»453 * 10“^ 

2,38 


Pilot plant 

data corresponds tor 



Catalyst life period 160" 

:'5hrs to. 1662.*Q 

hrs anti 


Pressure 26 kg/cm^ abs. 
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TABIE 4.3: CAICULAITIOH QB INTRINSIC CATALYST ACTIVITY 
A£ VARIOUS BEAOTI OH TEMP3HATUEBS 


No. 

Temperature 

°G 

Avg. value of 
experimentally 
determined rat 
constant ka^, 

hr~‘' 


1 

Oalci 
of ii 

HMa -VkO 4* ^ 

iiated value 
itrinsic 
c onstant 

-1 

sec 

® T 

' b 

“g jlcL u 

1. 

435 

0.3070 


105 

1.412 

* 

10-5 

9.06 

2. 

430 

0.2775 

* 

10^ 

1.422 

* 

10“^ 

7.715 

3. 

415 

0.2560 


10^ 

1.453 

* 

10“^ 

6.37 

4. 

412 

0, 2460 

•H* 

10^ 

1.460 

* 

10-5 

5.90 

5, 

410 

0. 2367 

* 

10^ 

1.464 


10-5 

5,50 

6. 

408 

0.2115 

* 

1o5 

1 . 468 


10“^ 

4. 47 

» 

7. 

405 

0,2160 

* 

10^ 

1.475 

* 

10-5 

4. 63 

8. 

400 

0.2242 

* 

1o5 

1. 486 


10“^ 

4.93 

9. 

395 

0.1948 

* 

10= 

1.497 


10"^ 

3.815 

lo- 

390 

0. 1840 


10^ 

1.508 

* 

10"^ 

3.42 

ll. 

385 

0,2080 

* 

105 

1,520 

* 

10"^ 

4. 25 

12. 

380 

0.1745 

* 

10^ 

1, 531 

* 

10-5 

3.09 

13. 

375 

0. 1790 


10^ 

1.543 

* 

10-5 

3.21 

14. 

365 

0. 1649 

* 

5 

10^ 

1.567 


10-5 

2.75 

15. 

355 

0. 1650 


10^ 

1.592 


10-5 

2.73 


Pilot plant data corresponds to; 

Catalyst life period 1114.0 to 1172. 5 hrs. 
and Pressure 21 kg/cm abs. 
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HI . B Effect of Pressure on Catalyst Activity 

Pressure serves to increase the activity of the 
shift catalyst. The increase in the activity wit^ pressure 
seems to he related to the availability of the internal 
surface area^^^\ 

To evaluate the dependence of the intrinsic rate 
constant on total pressure, an empirical equation has been 
fitted in the pressure range of 6-26 kg/cm^ abs. The same 
rate equation is assumed to be valid over the entire range. 

As all of the values of ka^ available for studying the effect 
of pressure are not exactly at the same temperature, a 
correction is necessary for the effect of varying temperature. 
Using the value of = 18,330 * 10^ gm cal/gm mole, the cal- 
culated values of kv^ have, therefore, been corrected for 
deviations of temperature from the level of 390®C. For fitting 
the empirical equation, the values of (kv^* P) have been 
calculated at 390°G, The results are shown in Table 4.4. 

A plot of log (kv^* p) against l/p is represented in 
Big. 4,3. A linear relationship is obtained which means that 
the relation between kv^ and P in the pressure range of 6-26 

O' , ■ ’ 

Kg/cm'^ abs, can be expressed by 

kViC< 2— j 

and the slope obtained from the plot gives 

K = 13.430 (when P = Kg/cm^ abs. ) 


(4.14) 



56 


I 


(Por extrapolation purposes, the same relationship may he 

2 

assumed to be Valid at least upto a pressure of 30 Kg/cm abs.) 


TABLE 4.4: EVALUATION OP THE PRESSUEB DEPENDENCE 
OP INTRINSIC CATALYST 


No. 

Pressure 

P,Kg/cm^ 

abs. 

Temperature 

°C 

Avg, value 
of experi- 
mentally 
determined 
rate constant 
ka ^ , hr“ ^ 

1/P 

Average 
Calculated 
value of 
(kv^*P) at 

390«C • 2 

Sqc"^ Kg/cm^ 

1. 

26.0 

380 

0. 2090*10^ 

0.0385 

110,82 

2. 

21.0 

380 

0.1970*10^ 

0.0476 

91.58 



390 

0.2001*10^ 



3. 


590 

0. 1930*10^ 




16.0 

390 

0. 1520*10^ 

0.0625 

62. 11 



590 

0.1610*10^ 



4. 


390 

0. 1592*10^ 




11.0 

388 

0. 1349*10^ 

0.0909 

48. 59 



380 

0. 1302*10^ 





375 

0. 1199*10^ 



5. 

6.0 

575 

0.6750*10^ 

0.1667 

18.44 


Pilot plant data corresponds to: 

Catalyst life period 266.0 to 578.5 hrs. 
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III.G Effect of Ageing on Catalyst Activity 

The activity of the freshly prepared shift catalyst 
undergoes a rapid decline during the first few hours of 
operation and thereafter, the rate of decline hecomes more 
gradual. The catalyst apparently sinters, thereby resulting 
in a decrease in surface area^^^\ This sintering phenomenon 
is the result of a combined effect of catalyst-age and reaction 
temperature^ But there appears to be no simple relation- 
ship between catalyst activity and these factors. It is 
rather difficult to individually assess these factors and 

correlate them with efficiency of catalyst ^ ^ 

f ) 

Mahapatra and coworhers' ' have computed the fall of 
catalyst activity with time for the H.T. shift catalyst using 
the data of a commercial reactor operating in the temperature 
range of 350-410®C. They obtained a straight line relation- 
ship on a logarithmic plot (shown in Appendix I?) of 
fractional activity against time of operation. Using this ■ 

I 

linear relationship, the dependence of the ageing factor 
(loss of original activity due to time on stream) on time 
can be expressed by 

Aj = (t >1) (4.15) 

where 

A^ - aging factor 
t = time on stream CdaV'^) 
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III.D Influence of H^S on Catalyst Activity 

The activity of a shift catalyst is lowered, by gases 
containing H2S which reacts with the iron oxide to form iron 
sulfide^^^\ In in<i^strial practice the CO - Conversion is 
performed without or with different concentrations of H2S 
as impurity, depending on the source of the feed gas. 

Although the effect of H2S on catalyst activity is 
quite important, not much work has been published in this 
field. To evaluate the influence of HgS, the use of an 
approximate sulfur- factor is necessitated in absence of any data 
for the piresent catalyst. Making use of the data from other 
sources^ a plot (shown in appendix IV) of sulfur factor, 
fg, against log (H2S) has been found to give a straight line, 
yhe variation of f^ with concentration of H2S can then be 
expressed by the relation 

fg • ^ 0,276 (H2S) + 1,127 (4.16) 

where 

(H2S) = concentration of HgS in ppm, (d^^y basis). 

At the point where fg=1, the corresponding value of (HgS) 
determined from the plot is 2.7 ppm which is somewhat uncertain, 
However, tli-is value may be considered as being practically 
equal to ?iero. 

Finally it may be mentioned that the use of eqn,( 4 ,l 6 ) 
for determining the sulfur factor is clearly an approximation. 

To calculate accurate factors would require experimental data 
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on the catalyst in the reactor, which are not available^ 

Combining eqns, (4.13) to (4.16), the final expression 
for the intrinsic reaction rate constant may be written as 


kv^ = Z* e 


■E./KT 


-F^/P 


• A^. f 
I s 


with 


Z^ = Z» 


-p^/p 


. A« f 
I s 


(4,17) 


Use of eqn, (4.4a) then gives 

-E /RT 

kvg = Z^. e ^ ^ — - . A^.fg (4,18) 


with 


Z 


0 



(4.18a) 


Evaluation of Z^ 

Use of eqn. (4. 13) and Eig, 4.1 has been pide to caiculnt® 


= 1,983 * 10^ sec”*'^ 

corresponding to the following conditions: 

t = (1605,5 + 1662, 0)/2 hrs, = 68 days 
P s: 26 kg/cm^ abs. 


yjj Q — 0, 4490 
(HgS) =s 60 ppm (dry basis) 


The values of yg q and (H2S) are available along with other 


pilot plant data, 


Using these conditions, 
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^ = 0.023 

= 0.716 

= 0.637 

o 

Then use of eqns. (4.13) and (4.17) gives 


Z' 


1.983 ^ 10° 
0.025 * 0.716 0.637 


1.895 * 10® sec" ^ 


The value of Z^ can then be evaluated from eqn. -(4. 18a), 


1.895 # 10 
0.4490 


,8 


= 4.22 * 10® sec"”* 


Eqns. (4. 17) and (4.18) can then be used to calculate the 
value of kv^ at a given set of conditions | yjj q, T, P, t and 

^ Vll • 

The value of E can be computed using eqns, (4,5) to 
(4.7) and (4.10) for a given catalyst size.t To give an idea 
about the magnitude of effectiveness factor, the values of 
E obtained at different positions in the catalyst bed are 
shown in Appendix IV), By inserting the value of E and kv., 
the apparent rate constant ka^ is calculated by means of 
eqns. (4.8) and (4.9). obtain th^ corresponding value of 
overall reaction rate, r, to be inserted in the material ", 
balance eqn.(4. l), use is to be made of eqn.(4. 4) with k^ = ka^. 
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lY, The Equili'briuin Oonversion 

The equilibrium conversion of GO can*"'''B^ evaluated 
from the value of equilibrixim constant for the shift reaction. 
The equilibrium constant, K, has been calculated from the 
thermodynamic considerations using entropy data and the 
relation thus developed is; 

In E = ^ 1^9998.220 + ( Aa .+ 10.414)T - T^- -4.^ 

- Aa.T. In (4.19) 

The details of the derivation of this equation are 
described in Appendix V, Some of the calculated values for K 
are listed in Table 4^5, These -are closely in agreement with 
the values reported by Lavrov'^ 

TABLE 4.5: EQUILIERIITM CONSTANT VALUES POR THE 


T^ »K 

fequilibrium Sonstan^ t 


M MIMi!Mi|i|Wi llllllilliiWIHiliM II II I IIWI ■'Ml I I I'M 

600 

28,410 

28.274 

700 

9.409 

9.4126 

800 

4.194 

4.2191 

900 

2.276 

2.2972 


V 


For computing the equilibrium conversion, the composition 
ratio Ey is related to the equilibrium constant by the 
relation^^^^i 
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K = K X Ky ' (4.20) 

Where 

An = change in the number of moles as a result of the 
reaction. 

K = composition ratio defined in terms of mole 

O' 

fractions 

yJoj % 

" yso ' %o 

Defining, 


^conv ^f , CO ~ ^5o 


K - illxgQz 
y " 


^conv^ ^conv^ 


^^f,G0 ~ ^conv^ ^^f ,'H^0~^conv^ 


K.- 


X 


' CO 2 ^P ^H2. 

( ^ ) h20 


(4.2la) 

(4.22) 


in which 

f = fugacity of the pure component at the temperature 
and pressure of the mixture. 

Por shift reaction, since there is no change in number of 
moles during the reaction. 

An =0 

Using the fugacity chart^^^^ ijo find'- the f/p ratio for each 
component, the values of iCc< have been calculated at several 
combinations of temperature and pressure in the normal operating 
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range of shift reaction- . The different ralues- of have 
been found to be very close to unity, so that 

K-. 1.0 

Hence eqn, (4. 20) becomes: 

™- K (4.20a) 

Eqns. (4.20a) and (4.21a) can then be combined together to 
give the value of and thus y^^ at any desired temperature 

from the given feed gas composition. 

V. Stoichiometric Relations 

The stoichiometric relations lead to the expressions 


for y. as follows: 

1 ■ . 

^f,aO' ^ (4*23) 

5^00 = yf,00 - yf.OO- ^ 

yoog = J'f.ooy yf,oo- 

^HgO = yf.HgO - yf,C0- (+-26) 

y^inert ” ^f, inert (4.27) 


Eqns, (4.1) to (4.10) and (4.15) to (4.27) thus 
constitute the simulation model for the shift reactor. 

The inlet temperature and flow rate of each component 
at the entrance of each bed (7=0) constitute the boundary 
conditions for the system, 

VI. Numerical Solution 

The simulation equations comprise a set of coupled and 
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nonlinear first-order ordinary differential equations vdiich 
cannot be solved analytically. The equations have therefore 
been solved numerically using fourth-order Runge-Kutta^^^^ 
integration technique starting at the bed inlet and proceeding 
stepwise, the computation cycle is continued until known value 
of the volume at the bed exit is reached. In all cases the 
maximum step size is chosen subject to the limitation of 
relative error specified (errors of less than 0 , 1 ^ in tempera- 
ture have been allowed). 

The solution to the simulation equations has been 
carried out on IBM 7044 computer using Fortran IV language 
for programming, 

VII ^ Comparison of the Simulation Model with Plant Bata 

A comparison of the model production with actual plant 
measurements is shown in Table 4,6, (A complete list of plant 
data has been presented in Appendix IV), Pig, 4,4 and 4,5 
show the calculated temperature and conversion profiles in 
the reactor beds. 

The agreement between plant data and model results is 
seen to be quite satisfactory. An exact matching of the two 
may not be possible in view of the errors in plant measurements 
and the experimental uncertainty in the kinetic data. 
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TABISJ 4»6: ^J^PAKSON OP REACTOR SIMULATION WITH 
PLAIT RElJlS^S ^ 


Parameter 

Simulation 


Plant Measurement 

1st bed 

2nd bed 

3rd bed 

1st bed 

2nd bed 

3rd bed 

^^composition 




» . 


of exit -gas 








62.56 

65.08 

65.51 

59.6 

63. 1 

64.6 

CO 

12.48 

4.90 

3,62 

14.2 

4.8 

3.6 

COg 

22.56 

27.78 

28.66 

23.0 

29.4 

30.0 

Ng+CH^^A. 

2.40 

2.24 

2.21 

3.2 

2.7 

1.8 

Temperatures 
oC (In/Out) 

'f 

351/475 

400/444 

406/414 

351/460 

400/459 

406/415 



2nd Bed jra 








CHAPTER 5 


OPTIMIZATION OF THE REACTOR 


The simulation model formulated so far is used in an 
optimization scReme to acRieve tRe maximum amount of reaction 
product from tRe existing multi-Red reactor, TRis is equi- 
valent to maximizing tRe final extent of tRe reaction witR 
a given feed state. For fixed amount of catalyst in eacR bed, 
tRis reduces to tRe determination of tRe inlet temperature, 
for eacR bed wRicR will give tRe maximum final conversion of 
CO. 


TRus tRe objective function, wRicR is to be maximized, 
is tRe overall conversion, of CO and tRis can be expressed by 
(Ref, Fig. 2.2): 

^ <^n,L - (5.1) 

n=1 


WRere 


X 

X 


N = number of beds in the reactor 
Q = degree of conversion at the inlet of bed n 
j = degree of conversion at the exit of bed n 


TRe optimum temperature profiles over the reactor are to be 
determined subject to the follovdng constraints, 
i) A fixed volume of each bed gives 





I ! 




n 


I. y 


f,CO 


"n,L 
n dx 


^rL,0 


50 


n = 1,2,. ..,N (5,2) 


,* 


\,L ^n,L 


( 5 , 3 ) 


This represents the natural constraint on the bed exit 
temperature which must lie below the equilibrium temperature 
corresponding to the exit composition, 

iii) In addition, the reactor temperature may be constrained 
below a maximum permissible value so that 

T ^ V (5.5a) 

Por the present reactor, the temperature may be limited to 
500®G^^^^ to prevent the deterioration of shift catalyst at 
higher temperatures. The temperature would exceed 500°0, if 
at all, only toward the end of the first bed. Therefore, 
the imposition of this temperature constraint may not be 
utilized in each and every case. 

Opt imi zati on-Method 

The solution of the optimization problem under consi- 
deration may be obtained either by a direct application of 

( 57 ’) 

dynamic programming' , or by using a direct search technique 
based on analytical criteria^ for the optimal policy. The 
dynamic programming approach may involve considerable compu- 
tation. The direct search approach, on the other hand, requires 
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numerical evaluation of the integral (Ref, Eqn,(2.6)) which 
involves the temperature derivative of the rate function. 

The experimental rate data must then he known with higher 
accuracy in order to find out this derivative. 

In the present work, another technique^ has there- 
fore been used which is based on a combination of analytic 
criteria for the optimum and a dynamic programming formulation. 
This technique is based essentially on the dynamic programming 
approach but utilizes a simple numerical search routine for 
single-bed optimization and the rate matching criterion of 
eqn. (2. 7) between two beds (n-l) and n* These search routines 
are described here briefly. 

Single-bed Outimization 

The optimal policy for a one-bed reactor with a fixed 
volume and given exit composition is determined using Modified 
Fibonacci search technique^^*^^. This amounts to a single 
variable search over reactor exit temperature within the 
specified range and finding the corresponding inlet conditions 
(by carrying out backward integration on the catalyst bed) in 
order to achieve the maximum conversion in the bed. The upper 
limit of the temperature range in which search is made is 
determined by the constraints represented by eqns, (5#3) and 
(5.3a) whereas a reasonable value for the lower limit is taken 
as 325°C, The equilibrium temperature to be used in eqn. (5.3) 
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may be calculated from eqns, (4. 19), (4. 20a) and (4.21) at 
the given composition. 


Rat e -Mat ching 

If the temperature ^ and stream composition q 
at the inlet of bed n are known, the rate-match criterion 
of eqn,(2,7) can be applied between beds (n-l) and n to 
determine the temperature, ^(^-1) L 

at the same composition (Ref. eqn, (2.7a)), This is accom- 
plished by using direct search over temperature. The procedure 
used for converging to the required temperature is as follows: 

A trial value for the temperature j, is selected 

which is higher than the temperature T^ ^ by some specified 
increment. The corresponding value of the reaction rate is 
calculated at the end of bed (n-1) and is compared with 
the value at the inlet of bed n. If the calculated value ' 
is higher, the process is continued using the same increment 
in temperature. On the other hand, if the calculated value 
is lower then the temperature has already passed through the 
desired point. At this stage the calculation is backed up by 
one trial and the process is repeated using a smaller incre- 
ment until the temperature is determined within desired 
accuracy. 


General Structure of the Algorithm 


Starting from the outlet end of the reactor, last bed 
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is optimized first using the single-bed optimization routine 
with an assumed value of overall conversion x,t t. Ibe 
remaining beds are subsequently optimized by means of the 
rate-matching criterion. If at the inlet of first bed, the 
relation 


■ 1,0 


= 0 


( 5 . 4 ) 


gets satisfied, the assumed value of Xjj is correct. Other- 
wise a new value is assumed (using a similar trial and error 
procedure as indicated for rate-matching) and the solution 
is repeated until eqn, (5.4) is satisfied within predetermined 
tolerances. 

While optimizing the reactor by this method, the upper 
temperature constraint of eqn. (5, 3a) has not been applied 
except for the last bed (for which single bed optimization 
method is used). As mentioned earlier in this chapter, the 
first bed exit temperature, ^ may then exceed the permi- 
ssible limit, 500°C), in some of the cases which is not 

desireable, The imposition of the temperature constraint 
in such cases becomes necessary to avoid excessive temperature 
and this thus eliminates the use of ratermatching criterion 
between first and second bed. 

In such a case, the problem is solved again (the 

actual of x„ ^ for this case must be very close to the one 
1T,L 

determined by the earlier method) by applying the earlier 
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method for optimization upto the second hed and then 
optimizing the first bed using single-bed optimization method 
which automatically takes care of the desired temperature 
constraint. An optimal ]Sf-stage policy is thus determined. 

The computations have been done on IBM 7044 computer 
using Portran IV language (the Program is presented in 
Appendix VI,). 


*** 


CHAPTER 6 


RESULTS MD liESCUSSION 

1, The optimum conversion has been calctilated for the 
existing feed composition and reactor load. Table 6,1 
presents the results of computation for this case. The 
actual plant data are also included for comparison. 

An increase of about yfo in conversion of CO from 
87.975^ to 91.145^ can be achieved provided the inlet tempera- 
ture to each stage is properly adjusted. These temperatures 
determined for the optimum reactor performance are, however, 
far from the existing ones, 

2, Por observing the effect of different variables on 
the performance of the reactor, several different sets with 
varying steam/gas ratio, dry gas feed rate and catalyst-age 
have been studied. In each case the overall conversion of 
CO has been optimized at the existing dry gas composition. 
Tables 6.2 to 6.4 show the values of these variables for the 
different sets and the resulting optimum conversions and 
temperatures. 

2.1 Variation in Steam/Gas Ratio 

a. Over the range studied, it can be observed from 
Table 6.2 that if the dry gas feed rate is maintained 
constant, the optimum CO conversion increases though the 
increase becomes less sharp with increasing stean/gas ratio. 
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This behaviour appears to be as a result pf a more favourable 
equilibrium for the increased ratio but a decreased contact 
time because the total flow rate becomes higher when more 
steam, is present. Thus beyond a certain limit, an increasing 
ratio is not much effective in improving the conversion as can 
be seen from Fig, 6.1, 

b. The adiabatic temperature rise is lower in case of 
higher steam/gas ratio. This can be attributed to the fact 
that the heat of reaction is utilized in providing sensible 
heat to a larger volume of gas-steam mixture, 

c. The sets of .optimum inlet temperatures show a 
consistent pattern. The temperature to each stage becomes 
higher when a higher steam/gas ratio is used. The reason for 
this behaviour may be that the heat of reaction raises the 
temperature of the larger volume of gas-steam mixture to a 
lesser extent, thus decreasing the catalyst activity. This 
unfavourable effect can be compensated by using a higher inlet 
temperature, 

2. 2 Variation in Dry G-as Feed Rate 

a. An increase in flow rate at a fixed steam/gas ratio 
leads to a decreased optimum CO-conversion (Table 6,3). This 
seems to be natural because a higher gas rate should retard 
the reaction due to a decrease in contact time. Moreover, it 
can be observed from Fig, 6,1 that the decrease in optimum 
conversion with dry gas feed rate 'is a linear one. 
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b. The optimum inlet temperature is higher at the higher 
dry gas feed ra,te. The higher inlet temperature vd.th an 
increased catalyst activity seems to compensate, in part, 
for the reduced degree of conversion, 

2.3 Variation in Oatalyst-Age 

The effect of time of operation beyond a period of 270 
days, has very little effect on either the optimum temperatures 
or on the conversion (Table 6,4), The optimum conversion is 
reduced by a negligible amount and the temperature variation 
is within 1® to 2°C only. However with passage of time, 
mechanical ’properties of catalyst pellets will deteriorate 
and ' necessitate replacement after about three years, 

3. Finally after studying the effect of individual plant 
variables, a two-level factorial design has been used to study 
the effect of interaction between steam/gas ratio and dry gas 
feed rate at the existing composition of dry gas. The combi- 
nations of the variables with the resulting optimum conditions 
are shown in Table 6,5. 

Experiments for this two level design have been carried 
out at conditions shown in Figure 6,2. The results are also 
included in the same figure. Following observations can be 
made. 

a. Increase in steam/gas ratio increases CO conversion. 

b. Increase in dry gas feed rate decreases CO conversion 
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c. The effect of interaction of above two variables 
is marginal because the effect of one is offset 
by the effect of the other. 

4. In almost all the cases studied in this investigation, 
it has been observed that the exit temperature from the first 
bed is always very close to maximum permissible temperature 
of 500°G. Moreover the optimum exit temperature from the last 
bed occurs within 36 to 52®C of the corresponding equilibrium 
temperature. 


■ABIE 6.5: EFEECT OE IWOIBIIACTION BETWEEN STEAI-l/GAS RATIO 
AND DRY GAS FEED RATE 
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Iniet steam/ gas ratio 



Response 

°/o CO -Conversion (X) 

Object 

Maximization 

Calculated 

values 

86-64 92-93 

93-12 

91-14 95-84 

Ratio effect 

C7as feed effect 

Ratio X Gas feed 
interaction effect 

(X 3 .X 4 -X 2 -X 5)/2 =5.50 

(X 3 *X 5 -X 2 -X 4)/2 =- 3.70 

(X 2 2 =0-80 


F ig, 6-2 - Calculated 'Effects' of the variables. 






CH&PTER 7 


COIGLUSIOI'TS ARP SUgGESPIOES 

A mathematical model has "been developed for a shift 
reactor which simulates the plant performance with reasonable 
accuracy. The results of optimization indicate that the 
conversion of the reactor studied can be improved by about 3^ 
by appropriately altering the temperatures, A further scope 
for improvement exists in principle by the use of a higher 
steam/ gas ratio. These changes must, however, be supported 
with the financial and operating aspects of the plant, parti- 
cularly, when the optimization of the entire shift conversion 
loop is desired. In that case the optimization has to be 
performed based on a plant profit function which includes heat 
recovery and steam separation costs together with capacity 
limitation on various units. 

The optimization-results obtained in this investigation 
will serve as a guide and a valuable starting point for 
optimizing the entire plant. 


*** 
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APPEHPIX I 

GOEDITIOIS POR OPPIxMIZATIOF OP MULPIBED 
ADIABATIC RBAOPOR 

The prohlen of finding minimum holding time to 
achieve a certain conversion is equivalent to that of deter- 
mining the maximum conversion obtainable with a fixed holding 
time in a given number of sections^ The two problems are 
duals of one another and the solution to both of these is the 
same. 


The holding time of bed n (Ref. Pig. 2.2) in a multi- 
stage sequence can be written as 

y. 


9 


n,I 

r- 


n 


y 




f 


(1.1) 


with 


here , 


n, 0 


r = r (y,T) 


9 = reactor holding time 
y = mole fraction of product. 

Por an N-stage system, the total holding time, 9^, becomes 

S’ y T 

— '^n,I 

/ f (1*2) 


N 

®T = "h?: 


n 


n, o 


1 1 yj 

Por the minimization of the total holding time for the 


sequence, it is possible to use analytical techniques to obtain 
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direct criteria for the optimal policy. The development 
shown below follows from the work done originally by 
The approach is illustrated for a two-stage sequence for the 
sake of convenience. The function to be minimized may be 
written as from eqn, (1.2), 


Q 


T 


J 


J 


1,1 

n 


1,0 


r 


y 


2,L 


2,0 


r 


( 1 . 3 ) 


Since there is no reaction in the interchanger, 
^ 1,1 ~ ^ 2,0 

y -j Q and y^ p are considered known, then 
®T = ®T ^1,1! '^Z,o'> 


The necessary conditions for 9^ to be minimized are found from 
the calculus. 




). 


^ TT a Ip 

■•’ 1,0 2,0 


■\ Cl 

^ "T'— 7p ) m ^ 

-x.O ,Ii 


= ( 


0 


A 0 


T 


), 


y P P 

1 , 0 ; 2,0 


(i.h 


On applying tnese ’conditions the following resultsi,,_ 

y. 


A ® 


T 




1 , 0 - 


1,L 




1 


A 


0 r 


A T 


dy = 0 


( 1 . 6 ) 
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Similarly, 


■d e 


T 


2,L 




■2,0J 

finally, 

t Cj — j 


^ 1 , L ;^ 1,0 ^ 2,0 


1 o r 
-.2 ’ 


. dy = 0 


(1.7) 


Jrp m 

^ 1 , 0 ;" 2,0 




1 


( 1 . 8 ) 


0 


*• ^1,1^ ~ ^2,0^ 


(1.9) 


The optimal policy must therefore sat'lrsfy the criteria re- 
presented by equations (1.6), (1.7) and (1.9). ^Phis approach 
may be extended to the case of an N-bed sequence with the 
following results; 


n,I 


^n^O 


P, T 


dy = 0 


n = 1,2, ,1T (1.10) 


and 


*ii,0 


) = r(y._o; 


n = 2, 5, 


wliich may be written as 
^(n-l),L “ ^n,0 


n = 2,3,.... ,N 


.,N 

( 1 . 11 ) 

(I.11A) 
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APPEMDIX II 

A. Effect of Pressure on Specific Heat 

Prom the specific heat data of Joseph^ for various 
gases, it has been found that the pressure dependency for 
gases CO, CH^^, IT 2 and A is negligible upto 40 atms. 

but in case of steam it is significant. The pressure-dependency 
has, therefore, been taken into account only in the case of 
steam. 

Table II. 1 lists the specific heat data^^^^ for 
steam in the pressure range of 1 to 40 atm. 


TABLE 

11 , 1 : 

SPEGIPIG HEAT OP 

STEAM 


Temperature 


Specific Heat, 

Cp/R 


°K 

1 atm. 

10 atm. 

20 atm. 

40 atm. 

550 

4. 535 

4.708 

5. 260 

6.967 

600 

4.391 

4.618 

4.921 

5.723 

700 

4.522 

4.627 

4.754 

5. 043 

800 

4. 665 

4.724 

4.792 

4.937 


C - G 

On a logarithmic plot of — 2n) against P or T, straight 

lines are obtained as shown in Pig. 11,1 and II. 2. The value 

(C -C ) is the difference of the specific heats in the real 
P Po 

state and the ideal state (Ps=1 atm.). Thus, for steam an 


equation of the form 
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(C - C ) = U • (II. 1) 

P Pq 

is valid and the following values have been evaluated for 
the constants U, Y and ¥. 

U = 0.63 * 

Y = -5.56 
¥ = 1.29 

The above equation is applicable at least in the temperature 
range of 600 to 800°K and for pressures upto 40 atm. 

Effect of Temperature and Pressure on Heat of Reaction 

The standard heat of reaction at 25°0 and 1 atm. pressure 
may be calculated from the heat of formation data^^^^ of 
Table II. 2 by means of the equation; 

pro ducts) “* reactants) (H*2) 


TABLE II. 2; STAN33ARD HEATS OP FORMATION 

(Conditions; 25°C, 1 atm, pressure, 
gaseous substances) 


Compound 

standard Heat of Pormation, g-cal/g.mole 

Ah° 

KjCs) 

0 

00 (g) 

-26,415.7 

oOgCg) 

-94,051.8 

0 

m 

-57,797.9 


Per the reaction 


CO (g) + H^O (g) > 002 (g) + ^2 
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the standard heat of reaction is calculated by means of 
eqn. (II. 2): 

= (-94,051.8) + 0 - (-26,415.7) - (-57,797.9) 

= - 9,838 cal/gmmole 

Taking 600°K and 1 atm. pressure as the reference state (this 

corresponds to lower range of normal operating conditions 

for H. T. shift reaction), the following equation may he 

applied for determining the effect of temperature on the Heat 

of Reaction at constant pressure^ ; 

T 

AHj^(T,P) == AHj^(T^,P) + j ACp dT (11.3) 

T 

0 

where AC = change in heat capacity of entire system at 

constant pressure. 

^ *^p^ pro ducts “ reactants 


The heat of reaction at 600°K and 1 atm pressure is 

n T^600 

AHj^( 600°K, 1 atm.) = A + P Ac dT 


Tq=298 


(II.5a) 


kt atmospheric pressure, the gas is ideal and the molal heat 
capacities may he represented hy, an empirical equation of the 
form 

C a + hT + cT^ (II. 4) 

Then, 

Ao =Aa + AhT + AcT^ 

P 

Using eqn. (II, 3a), the integration yields 


(11,5) 
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^Hj^( 600°K;, 1 atm.) = A H° + [Aa(5}-T^) + (T^-T^) 

+ -§^(1^-1^)] (II. 3b) 

in which = 298°K 

0 

T = 600°K 

~ -9838 cal/gm mole. 

Using the values of empirical constants from Table 4.1, 


A a 

= '"CO 2 

•f 

■" ^00 ‘ ^H20 

= - 0.201 

Ab 

- 



= 5.493 * 10 "^ 

Ac 

= °C02 

+ 

" °G0 " 

= -2.7177 * 10 


Substituting these va,lues in eqn. (ll,3b), the heat of reaction 
at the reference state (600°K, 1 atm.) is calculated to be 

(600°K, 1 atm.) = - 9328 Oal/gmmole, 

To evaluate the effect of pressure on heat of reaction at a 
constant temperature of 600®K, the principle that the enthalpy 
is a point function results in the equation^ ; • 

AHp_(T,P) = AHj^(T, 1 atm) + Chj,(T,P) - Hp(T, 1 atm)] 

(II. 6) 

v/here H is the enthalpy and subscripts R and P refer to 
reactants and products, respectively. 

Use of Eqn, (11,6) for the shift reaction then gives 

AHj^(T,P) = AHj^(T, 1 atm) + |[qq^(T,P) + Hjj^(T,P) - Hqq^(T, latm) 

- (T,latm)J 

-(HQo(Itl) + Hjj q(T,P) - HQQ(T,latm)-Hjj^Q(T,1atm3 

(II. 6a) 
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Por tlie gaseous reactants and products, tile enthalpy change 
in shifting from the lower to the higher pressure is evaluated 
using the data^^"^^ given in Table II. 5. 


TABLE 11,3: EICTHALPY OF BEAGTMTS MD PRODUCTS 


Pressure 

atm. 

Enthalpy at 600 °K, 



00 

H2OU) 



1 

7.7546 

8.965 

9.802 

7. 6084 

10 

7.7528 

8.811 

9.764 

7.6147 

40 

7.7486 

8.216 

9.635 

7.6359 


= The aero reference point of the enthalpy, taken 
as the internal energy of the ideal gas at 

absolute zero. 


R = 1.987 Gal/ijiole "K 

T^ = 273.16°K 

A plot of- AH j^( 600°K, P atm.) against P(in atm) on a semilog 
graph gives a straight line. The relation obtained for the 
variation of vdth pressure at 600°K is 

AHj^( 600°K, P atm) = - (0.933 * 10^^) exp (-0.9354 * 10"^P) t 

(II. 7) 

To evaluate the heat of reaction at any temperature and 
pressure in the normal operating range (600° - 800°K, 1-40 atm.), 
Eqn.(II.3) may be used with T^ = 600°K. The value of 
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(600°K, P atm.) may be substitute! from eqn, (11,7), The value 
of above atmospheric pressure may be determined from the 

follovdng eqn, ^ 

AG^ = Aa + Ab.T + Ac.T^ - UP^ (11.8) 


TABIiE II. 4: CALCUL;i.TED VARIATION OP ENTHALPY 

DBPARTUEE AND HEAf"W”EBr<3TOrWITH 
PRPiSSURE ( TtoERATM 600 


Pressure, 

P atm. 

Enthalpy Departure at 6O0®K,cal/mole 
AH(600°K,P atm. )=H(600 oK,P atm)- 
H(600«K, 1 atm) 

■SeaF'pT 

Reaction at 
600«K and 

P atm. Pres su] 

CO 

HjOCg) 

COg Hg 

" AHp(6Q0®IC, 
atm) 
cal/gaanole 

1 

O 

0 

0 0 

•9328 

1 

-0.978 

-83. 600 

-20.64 3.424 

-9260,64 

40 

-3.260 

-407.000 

-90.700 14.94 

-8993.5 


The final expression obtained after evaluating the 
integral term in eqn, (II. 3) is given by 


AHj^(T,P) = AHj^(600'»K,P) +[Aa.T + + -^4 T^ 


■with 


= 672.47 

Cg = -13.82 * 

= - 2.96 


- C2p''t^+^]4.C3 


P 



1 

Ui.9) 
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APPENDIX III 

IBRIVATION OP 'IHB PSEUPO-PIRSO? ORDER RATE EXPRESSION 


Por determining the reaction rate of CO conversion, 
the second order expression may he used for the reversible 
reaction 

^2 

CO + H O (g) - --- ' ' COp + 

^2 

The overall rate of reaction can be expressed as: 


r 


= *=2 ^00 




where ; 



(m.1) 


k^, k^ = second order rate constants for the fgrward 
and the backward reaction respectively. 

Since catalysis does not change the equilibrium but only 
increases the reaction rate in both directions, the relation, 

K s = ^®^ween the rate constants and the equilibrium 

constant is valid also for the heterogeneous reaction. 

Equilibrium constant K is related to the equilibrium 
composition by: 


K 



X y|^ 


(III, 2) 


Substituting kk = kp/K in eqn.(III,l), we have 

^COp 

r = kj [yoo- yHjO K J 

Combining with eqn, (III,2) 


(III. 3) 
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= '"2[5^00-yH20- 


CO, 




nd' W' 


Since H20(g) is present in large excess, its concen- 
tration may be ass-umed to remain constant and equal to the 
mean value y^ q. 

The eqn, (III-, 4) then becomes 


with 


= ^^2 %0 &C0 - " ygo] 

^002 ySjO 

” ° ySog ■ yfj ■ 5^ ^ 


(III. 5) 


The more the reaction approached equilibrium, the more 
b approaches unity, and it would become equal to one only 
for the equilibrium state. The equilibrium concentration, y^^, 
for this reaction is only a small fraction of the original CO 
concentration in the inlet raw gas-steam mixture so that the 
error of taking b=1 would be negligible in the initial as 
well as in the final stage of the conversion. 

Putting b=l, the final expression for the rate of 
reaction then becomes: 

r =: k,j y^Q - y^Q ) (III, 6) 

with k^— ^2 0 

Equation (III, 6) is a pseudo first order rate equation in CO 
and k,j may be termed the *pseudo first order rate constant” 
which depends, naturally, on the actual concentration of 
HgOCg), 



APPENDIX tv 


MISGELLAEBOUS RELEVANT DATA 


TABEE IV. 1 : SULFUR EACTOR DATA 


( 35 ) 


Sdncen^ratibn’ of in the gas 
(dry Basis) 


Sulfur factor, f^ 


1000 


1.000 

0.650 

0.575 

0.525 

0.490 

0,460 

0.440 

0.419 

0.400 

0.590 

0.580 

0.299 

0.259 


2000 




«! 
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TABIB IV.2: CATALYST EBEECTIVENESS BACTOR 
GALCULAT]^ ALON(^ THE lItIi&TII Of 



THE REACTOR 

BED* ■ 

Srac-tional 
bed len^cth 

Temperature, 

Effectiveness factor 

•0 

624. 1-6 

0.2673 

0.1 

637.34 

0.2364 

0.2 

651. 13' 

0. 2085 

0.3“ 

665.35 

o 

00 

o 

0.4 • 

679.65 

0.1630 

0,5 ■ 

693.66 

0.1453 

0. 6 

706.97 ' 

0.1308 

0.7 

719.27 

0.1190 

0.8 

730.30 

0. 1097 

0,9 

739.95 

0.1023 

1,0 

748, 18 

0.0965 

*1st 

bed with the feed conditions as given 
in Table IV. 3. 

TABLE 

IV. 3: PLANT DATA 

FOR THE SHIFT REACTOR " 


I* Catalyst 

Types ECI H,T. shift catalyst 

Catalyst size and shape ; 6 xm ^ S nun , cylindrical pellets 
Ai^e of the catalyst: 270 days 
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II, Reactor 

Ro. of beds; 3 

Reactor bed dia, ; 9’6» (I.R. ) 

Catalyst distribution: 

1st bed; 154 cuft. 
2nd bed: 308 *’ 

3rd bed: 308 '» 



3 rat ure 


(ih/Out ) 


1/460 


400/459 


406/415 
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APPENDIX V 

EQUILIBRIUM CONSTANT POR WATER GAS SHIPP REACTION 

The relationship between the standard free-energy 
change A P® and the equilibrium constant K is given by the 
expression: 

AP° = - RT In K (V. 1 ) 

The value of Ap° can be determined from the expression 
Ap° = AH° - T a S° (Y, 2) 

where 

AH° = Standard heat of reaction 

= Standard entropy change for the reaction. 

The value of AH*^ is evaluated from the data for standard 
heats of formation (Table 11,2): 

AH° (298°K) = (-94052 + O) - (-26416 - 57798 ) 

= - 9858 cal/gmmole 

and the As° Is evaluated from the Entropy data of Lavrov^ 
given in Table V. 1. 


TABLE V. 1 : STANDARD VALUES OP ENTROPY 


(Reference temperature » 25 

Compound 

S° cal/mole =>0 

OOg(s) 

51.069 

H, (e) 

c 

31.2110 

00 (g) 

47.216 

HjOCg) 

45.106 




85 


Thus 

As° (298°K) = (51.069 + 51.211) - (47.216 + 45.106) 

= -10.042 cal/(gmmole)( °C) 

The temperature dependence of equilibrium constant can 
be determined using the data for heat capacities (Table 4.1) 
of reacting substances with the following equations: 

tO / m \ A ttO I 


T 


and 


A h“(T) = AH°(Tq) + r AC dT 

T- ^ 
o 


T 




As°(t) = As°(t^) +j ru 


(II. 3a) 


(V.3) 


T 


The quantity Ac^ is given by equation (II. 5) and from this 

T 

f ACp dT = Aa.T + ^ ^ ^ _ -,60.220 

t'^=298°K 

rn A C 

i p A 

dT = Aa, In T + Ab.T + . T^ - 0.372 

T'^=298°K 

Substituting these values and combining eqns,(V.2), (V. 3) and 
(II. 3a), the fiml expression obtained for Af° is 


Ali'°(T) - - 9998.220 + (Aa + 10.414)T 


Ab 


T' 


c 


. T^ - Aa.T. In f 


2 ' b 

an d i n n e rt i rig this in e qn « ( V , 1 ) give s 


(V.4.) 


In K (' zj ij r"' ) t^998,220 + (Aa + 10.414) T 

. ^ 0}2 ^ ,j3^ Aa.T. In T^ (7.5) 



APPESTDTy -tr-]- 


COMPUTEr PROGiUWMB 
for 


OPTIMIZATION OP REACTOR 
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C* »*•»■* ** Jf-*- **# 

., C 

OPTIMIZATION-PROGRAM 

C 

C»- «■* **«■ » »■«■■»■«■«•■!!••«•***«•*# 

c 

C«»» SIMULATIO!>l DPTIMISATION OF WATER GAS SHIFT REACTOR 

C*»* UPTIMISATIOM TO GET MAX«D'-GR^i5 OF CONVERSION IN MULTIBcD ADIABATIC 
C*** REACTOReTHE NO«QF BEDS AMO VOLUME OF EACH BFD IS SPECIFIEDo 
C«**- mater IAL-FMF;P.GY BALANCE ECMSo (TWC FIRST ORDER SIMULTANEOUS 
C**# DIFFERENTIAL EQNSo ) APE SOLVED USING RUNGE-KUTTA METHOD® 

C«** RATE EQN. IS PSEUDO FIRST ORDER CBTAINED BY SIMPLIFYING 
C*** SECOND ORDER RATE EKPRE SS ION ( LAUPI CHLER , MOE ) 

; C«** FOR FCI H»T«SHIFT CATALYSTtCHAMDRAS PILOT PLANT DATA)® 

^ C“»- NOTATION FDR I'JPUT DATA 

C~»« PERC{NC) = INLET GAS COMPOS IT I ON { VCLUMF PERCENT DRY BASIS) 

C-*- COMPONENT NO® AND NAME IS AS GIVEN — 

C-*~ 1=H2, 2=CQ, 3=Ca2T 4=CH-^, 5 = N2, 6 = A, 7=H20 

FDGSn.)=INLET DRY GAS FLOW , RSG { 1 ) =STEAM/G AS RATIO INLET 

C-*- P'=PR»:SSURFt DP=CATALYST SIZE 

C*** UNITS— TIMDAY=TIMF IN OAYSCON STREAM )-—( GREAT ER THAN 1 DAY) 
fC«** UNITS— H2SPPM=CONCtUOF H2S IN PPM^(DRY BASIS) 
i C-»*» P = KG/CM2 ABS 
lC**» DP=PARTICLE DIA=CM. 

[C ^ ##*### # #### ### *## ■#•#####. 

I DIMENSION PERCdO ) ,FC(aL.v)»CQNV{ 1 ),FDGS(25) tRSG(25)tFTS(25) 

I OIMEMSIOM YCO'lVd. ) ,YF( li )tY(l-y'' ).TCS(25) ,X{25) , VOLUME!?. ) 

I DIMENSION KEOPTMtl'OtTIOPTMd"') tTEOPTMd?) tTEEQLMd.)) 

' DIMENSION KEd* d.TEdO tXId. )»TI(iO) ,RATI(lD),XIOPTM(l;:«) 

DIMENSION ACP( 7),BCP{ 73,CCP{ 7) 

TKO=273ol6 
NC0MP=7 
NS=1 
NCl = l. 

NC2=NC1+1 

NC6=MCl+5 

NC7=NCl+6 

data ACP/ 6» 945 t 6« B50 j 6o 3? 9» 3o 204 t6a 45Tf 4 o 97 » 7c 1 36/ 

DATA BCP/-3»196E--3» IcSllE-Btl o 14E'»3 1 1S«41E”*S ,1 *3 B9E“3 » » 2c 64E-3/ 
DATA CCP/t...4757E“6t »267Sfc-6,-3c4*5C-6,-4.48E”6t- '*» :‘69E-6t'.'c...» 
10cL459E-6/ , 

READ 146, P, DP,N3E0 
146 FORMAT! 2F11#2, 111) 

READ 143, !PERC(NC)tNC=NC1tNC6) 

143 F0RMAT(6Fll«2) 

READ 144, FDGS(NS) ,RSG(NS) 

144 FORMAT! 2Fllc2) 

READ i4S, (YOLUM^tMB ) ,NB=l,N3eO) 

145 FORMAT! 3Fllc3) 

RFAD 147, H2SPPM,TIMDAY 


87 


147. F0RMAT(2F1I»1 ) 

, READ 148, TI MINT, NTI MEL 

148 FORMAT! Fll, 1,111) 

PRINT 445 

PRINT 455 
PRINT 999 

999 F0RMAT{/4X,»riPTIMIZATI0f! OF TROMBAY SHIFT REACTOR ““-AT DIFFIRIMT 
ICOMSINATIDNS OF STEAM/GAS RATIO AND DRY GAS FLOW*/) 

PRINT 455 

C cx, * E. * — ^ „ jj. ^ ^ ^ s,|. 

PRINT 15 ^) 

15'- FORMAT! //AX, ♦INPUT DATA») 

PRINT 156, ! PEPC!NC) ,>1C = NC1,NC6) 

156 FORMAT!4X,*PLRC!NC)*,5X,6F13,2) 

PRINT I57,\IS, FDGS!NS) ,RSG!NS) 

157 FORMAT!4X,*NS=«-, 12, 3X,*FDGS( NS) = *,Fl2o2,5X,*RSG!NS )=*,F6.2) 

PRINT 161,! VOLUME !NB) ,M6=i,N8ED) 

161 F0RMAT!4X,»VDLUME!N3)«-,5X,IXF",4) 

PRINT 165,7,DP,N8ED 

165 F0RMATI4X, *•?=•», F6®7,5X,*DP=*,F5c2,3X,»NBFD=*, 12) 

PRINT 166,H2SPPM,-TIMDAY 

166 F0RMAT!4K,#H23PPM=»,F6«l,5X,*TIMCAY=»,F7.l) 

► PRINT 167,TIMIMT,'NTIMEL 

167 F0RMAT!4X,»TIMINT=»,F7o l,5X,*NTIMEL=*,I3) 

PRINT 17!^ 

17* FORMAT !4X, ♦END OF INPUT DA,TA») 

DO 600 NR SG AS =1,3 
DO 1«:< MC=NC1,NC6 

15' FCINC)=FDGS!NS)*PERC{NC)/IC-D«0 
FC!NC7)=FDGS!NS)«RSG(NS) 

FTS!NS)=FD3S(NS )+FC !NC7) 

FTSR=FTS!NS) 

F=FTSR 

DO 15 NC=NC1,NC7 
Y! NC)=FC! NC )/FTSR 
15 YFINC)=Y! NC) 

C*** FEED COMPOSITION CALCULATIONS ARE OVER. 

C # # # # # # » . # # # ' # ^ «• # '# # * * * * * # » # # 

DO 560 NTIME=l,NTIMFL 
PRINT 456,TIMDAY 
PRINT 5„F 

5 F0RMAT!//4X,*V0LUMETRIC FEED RATE OF WET GAS=F=*, F13. 2 ) 

PRINT 65 

65 F0RMAT!///4X,«-C0NVERT0R INLET GAS COMPOSITIONIFROM INPUT DATA)*) 

GO TO 79 

19 XF=f:'.f’ 

VI=m.-n 
PRINT 445 
PR FIT 455 
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C«** START OF OPTIMISATION CALCULATIONS^ 

SINCE XP(DEGR;-E OF CONV„OF CO IN PRODUCT GAS) IS UNKNOWN UijTIL 3PT I MI S/.T ION 

CALCULATIONS ARE COHRLFTRo SO START WITH AN ASSUMED VALUE OF K P, CALCUL4T c 

— - CORRESPONDING VALUE OF XFCFOR FSED)AND COMPARE THIS WITH GIVEN XFo 

— TRIAL AND ERROR METHOD IS REPEATED TO OBTAIN XP UNTIL 
CALCULATED XF=GIVE:N i'.Fo 

**»*«■«*»(■**»» 

— ** CASE 1-“CNSTRN='; « (IN THIS CASE » CONSTRAI NT OF MAX«TEMP«OF j , D = GR€n C^ 

IS NOT APPLIrD«MULTI3ED REACTOR IS OPTIMISED USING RATH-MATCH CRITERIA)c 

CASE 2— CNSTRN = 1« (IN THIS C ASE t CONSTRA I NT OF MAXoTFMP^OF S'" DEGREF Cc 

— - IS APPLIED FOR 1ST BrDoSO RATE-MATCHING CRITERIA IS MOT USED FOR 

— - OPTIMISING 1ST BED OF MULTIBED REACTOR^IST BED IS OPTIMISED US I (G 
— “ SINGLE BED OPTIMIZATION ROUTINE), 

*'«•***»» »«•*** «■ *#■«•** 

TMXLMT=5T<’of +273o- 
CNSTPN=,'.. 

XESTRT='’'. 96' : 

OFLTAX = ’'-o'"l.'v 
KX='' 

NB=N8ED 

XE ( MB ) =XE STRT-DELTAX 
A 15 XEN=Xr-(NB) 

XP=;<F(NB) 

VS=VOLUME (NB) 

KX=KX+I 

25). CALL SBOPTM(KHN,THN,P,YFt ACP,BCP,CCP,F,XPtDP,XIN,TIN,VItVE, 
iTIMDAY»H2SPPM) 

TE(NB)=TEN 
XI(NB)=XIN 
TI{NB)=TIN 

IF (NB ED-1 )3)'''' »3f 0,475 
A75 IF(CMSTRN )43",26i;,A8 
ASt.'^ IF {NB-2)3rT' 1485,260 

A83 XE(NB-1)=XI (N3) 

NB=NB-1 
XEN=XE(NB) 

GO TO 25D 

260 CALL RATMCH(XIN,TIN,P,YF,XP,DP,RATIN,TENl,TIMDAY,H2SPPM) 

RATI(NB)=RATIN 
TE(NB-l)=TFNl 
XE(NB-1)=XI (NB) 

NB=NB-1 
TEN=TE(NB ) 

XeN=XE(NB) 

V&=VDLUME (NB) 

CALL MATENG(XEN,TEN,P,YF,ACP,3CP ,CCP,F,XP,DP,XIN,TIN,VI fVfc, 

1TIMDAY,H2SPPM) 

TI (MB)=TIN 
XI{NB)=XIN 
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IF(N8-1)3''" ,3' r ,47? 

SHf' DO 325 MB=l,NeCD 
XIOPTM(HB }=XI (N3) 

XE0PTM(N8 )=X=(NB) 

TIOPTM(NB )=TI(M8) 

TEOPTM(NB ) = TI-{N8) 

323 CONTINUE 

PRINT 345,S13r:D,Xe(NBE0) 

345 FORMAT! /5X, ^ACSUMifO VALUE OF ^EGREE OF CONVERSION AT REACTOR EXIT 
l=Xc(^^,I2,»)=*,F3o5) 

PRINT 3‘4'fc 

346 FORMAT ( /TX,*NB*, 9X, *XIQPTM ( N8) *, 9X,»XFOPTM ( NB ) » ,9X, *TI OPTH( NB ) *, 
:9X»*TE0PTM( NB)*) 

PRINT 347, ( NB,XIOPTM!NB ),;.;E0PTM{ NB ) , TIOPTM! NB ) , TEOPTM { NB ) » 
1NB=1,NBED ) 

34? F0RMAT{6X,I3,2F19c4,2F19o 2) 

42f/ XI1DIF=XI0PTM(1 )“XF 
XIEROR=ABS( XIIDIF) 

IF(XIEROR-C<.; I' )45 T425t425 
425 IFtXIlDIF )43. ,450,43" 

430 OELTAX = DELT «250 
IF{KX-l)41it,^I' »A40 
435 MB=NBED 

XEOLD=XE(NB) 

44* , NB=NBF.O 

XE(NB)=XE0LD-DELTAX 
GO TO 415 
45' PRINT 455 

455 FORMAT! /IXt 12/ ( IH*) /) 

DO 33'' NB=1, ,NBED 
XEXITB=XE(MB) 

TEEQLMINB )=TEQ! YFjXEXiTB) 

33i: CONTINUE 
PRINT 445 
445 FORMAT! IHl) 

PRINT 455 
PRINT 46D 

46!: F0RMAT!/5X, 2H**,*FINAL OPTIMUM CONDITIONS FOR REACTOR BEDS* , 2H»*/ ) 
PRINT 456,TIMDAY 

456 FORMAT! 7X,»TIME ON STREAMilN DAYS) =* , F7»l / ) 

PRINT 223tF0GS! 1 ) 

22* F0RMAT!4X,»INLET DRY GAS FLOW=» , F12,2 ) 

PRINT 225tRSG! 1) 

225 FaRMAT!4X ,»INLET STEAM/GAS RATIO=* ,F6« 2 ) 

IFCCMSTRN )51f:,50P,5Xn 
5?.:i'* PRINT 5u5 

5; 5 FORMAT! 7X fifCONSTRAINT OF MAXoTEMP«OF SwODEGRcE C« I S NOT APPLIED®*} 
PRINT 506 

5*’6 FORMAT! 7X,*MULTIB5D REACTOR IS OPTIMISED USING RATE-MATCHING CRITE 
IRIA*/) 
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GO TO 520 

51‘-' PR INT 515 

515 FORMAT (7K,*C0NSTR&I NT OF MAX^TEMPoOF 5L",> DEGRHE Ce IS APPLII-D FOR 
ilST BFD*) 

PRINT 5i6 

516 F0RMAT(7X, »RATE»MATCH ING CRITERIA IS NOT USED FOR OPTIMISING 1ST 
•-BED*/ ) 

52f CONTINUE 
PR INT 348 

34,3 rQRMAT(/6X,*BED»»llX,*INL5T»,15Xt*EXlT*,16X»*lNLFT*t 
114X,^^!-XIT«-, :: X,«EQUILIBRIUM*) 

PRINT 349 

349 F0RMAT(6X,»Mn«,*,9X»«-C0NVERSI ONfr , 9):,*C0NVERSI 0M» tl3X ,»TF;MP» » , 
114X»*TEMP.» ,1 X,*r:XIT TEMPo*) 

PRINT 35D 

3 5' - FORMAT ( / 7X, »N8* , 9X, *X lOPTMt NB ) *»9X,*XE0PTM ( NB ) » ,9K , *T lOPT M( MB)*t 

19X,*TE0PTM( NB)#»9X,^^T’--EQLM(N5)*) 

PR INT 351 » ( MB,XIQPTM( NB ) t XEOPTM < NB ) ,TI0PTM ( N8) t TEOPTMC N8) t 
ITESQLMINB ), MB = l,NBSO) 

3 51 FORMAT { 6X , I 3 , PF i9o 4 , 3F 1 9o 2 ) 

C •»«•■»*’»***•»■ *»*■»«■•»»• «■ 'it **■■*»•»«•»■ **»«■»«■«■#«■»**•*# ■»■#*»# 

C*** FINDING COMPOSITIONS FOR DIFFERENT BEDS 
NB = 1 
NS = 1 
NC1 = 1 
NC2=NC1+1 
NC6=MCi+5 
NC7=NCl+6 

25 X(NS)=XIOPTM{NB) 

TC S J NS ) =T I D PTM ( '^1 B ) - TK Q 
X( NS + l)=XE3PTM(;iB) 

TCSINS+l) =TEOPTM CNB )”TKO 
XEXIT=X(NS+1) 

GO TO 1?» 

It'C NCl=NCl+7 
NC2=NC2+T 
NC6=NC6+7 
NC7=*NC7+7 
NCLi=NCl-7 

lt'2 CALL MASS(NC1,YF,«XSXIT,.Y) 

YCONV(NB) =YF(2)*tX( NS+1)-X(NS) ) 

33 C0NV(NB)=YC0NV(N8)*FTSR 
DO 35 NC=NC1,NC7 

35 FC (NC )=Y( NC)*FTSR 

FDGS(NS‘+1 ) = FTSR“PC{NC7) 

DO 37 NC=NC1»NC6 

37 PHRC(NC)=1DC,G*FC(NC)/FDGS(NS+1) 

PRINT 7(4 NBtCONVtMR ) ,NS 

7- FORMAT( // /4X» IIHBED NO® =MB= 1 2 1 16X, 2t}HCONVERSION=CONV ( NB ) = F9« 2 1 
15X,*STRGAM !‘iru=NS=*tI3) 
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PRINT 75,TCS(NS) ,'TCS(NS+1) ,X{NS) ,X{NS + 1) 

75 FORMAT {4X t-l 6 HTPMR 0 1 N=TC S ( NS ) =F7« 2, 6X,20HTEMP« EXIT=TCS ( NS+1 ) =F9, 2 1 
15X,»X{NS)=*, F7»4,'FX,*X{NS + 1) =*,F7o4) 

PRINT 95, NB,VOLUMF{ MB) 

95 format { /4X,*T0TA.L VOLUME OF BED* ,1 2,»=* , F7«,3 ) 

PRINT 96, MB 

96. FQRMAT(/3X, lH*,*GOMPOSITION OF GAS EXIT BED*,I2,1H*) 

79 PRINT 80 

8 i' FORHATt /4X,9HC0MRnNENT,.9X,4HFL0W,l XtllHMOLE FRACN« , 6 X, 13HPERCENT 
ICOMPo ) 

PRINT 82 

8 2 FORMAT ( 7X ,2 HNC , 1 1 X, 6HFC { MC ) , 1 IX , fiH Y { NC ) , IIX , 8HPERC CMC)/) 

84 PRINT 85, (NC, FC(NC) , Y ( NC) , PERC ( NC ) ,NC=NC 1 ,NC 6 ) 

8 5 F0RMAT{6X,.I3,irx,'F9o2,9X,F7*4,12X,F6*2) 

NC=NC7 

PRINT 88 , NC,FC(NC), Y(NC) 

88 F0RMAT{6X,I3,7X,.FI2o7,9X,F7,4> 

89 IF(NCl-l)i9,19,90 

9 C I F ( N B“ N BE D ) 9 2 , 6 0 6 f.' 

92 NB=NB+1 

NS=NS+2 
GO TO 25 
CONTINUE 
PRINT 455 

IF(TFOPTM (1 )-*TMX.yMT)55S ,550,525 
5 2 5 CN 5TRN = 1 « 

XESTRT = X£0PTM(N8ED) +-’o015f 
0ELTAX=C.r'>25 
PRINT 445 
PRINT 455 
GO TO 410 
55C CONTINUE 
PRINT 445 
NB = 1 
NS = l 
NCl = l 
NC2=NCl+l 
NC6=NCl+5 
NC7=NCl+6 

TIMDAY=TIMDAY+TrMINT 
56t. CONTINUE 

IF(NRSGAS-2)625,650,6C0 
625 RSG(NS)=1»5C 

FDGS(NS) = 22 5f'(,)*CC 
GO TO sm 
65C RSG(NS)=1«50 

FDGS( NS) = 13500.00 
6 "- CONTINUE 
STOP 
END 
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SUBROUTI^JE SBQPTM(XE»TEOPTHt P,YF »ACPtBCP,CCP,F,XP,DP,XIOPTM» 

ITIOPTMtVI rVFtTIMCAY,H2SPPM) 

— THE PURPOSE OF SUBROUTINE SBOPTM IS TO OPTIMISE A SINGLE BED OF THE 
REACTOR.VALUE OF BED VOLUME AND XE (DEGREE OF CONVoOF CO AT BED EXIT) 

— IS GIVEN, THIS SUBROUTINE FINDS TEOPTM ( OPTM. EXIT BED TEMP, ) t 
— - TIOPTMIOPTM, INLET TEMP) ,XIOPTM( OPTM. INLET DEGREE OF COMVoOF CO) 

TO GFT MAX, CONVERSION IN THE BED. 

TEMP, T IS IN DEGREE K, 

MAXIMIZATIOil USING MODIFIED FIBONACCI SEARCH 
DIMENSION ACP( 7),BCP( 7),CCP( 7),YF(iO) 

TEQe=TF.Q{ YFtXE) 

TMXLMT = 5ri;\ 1 +273,' 

K = 1 

TELOW=325,D+273,0 
IF(TMXLMT”Ti"QE) 5,'6,6 

5 TEHIGHsTMXLMT 
GO TO 7 

6 TEFIGH=TFQE 

7 CONTINUE 

1>' TrRAMG=TEHIGH-»TEyOW 
TELOCT=Ci, 38 2T*TERAN G 
TE1=TELOW+TELOCT 

C*** AS WE ARE PROCEEDING FROM LAST STAGE TO FIRST STAGEtEXIT CONDITIONS OF 
^C«** REACTOR IN THIS ARGUMENT LIST CORR-SPOND TO INLET CONDNS, IN THE 
C*** ARGUMENT LIST OF SUBROUTINE MATENG, 

CALL MATENG(XE»TE1iP,YF,ACP,BCP,CCP,F,XP,DP,XI1,TI1,VI,VE, 
1TIM0AY,H2SPPM) 

XCONVl=XE-KII 

TE2=TEHIGH-TELOCT 

CALL MATFNG{XE,T82*P>YF,ACP,.BCPtGCP,F,XP,DP,XI2,TI2,VI,VSf 
1TIMDAY,H2SPPM) 

XCONV2=Xe-*KI2 
2U TEDIFF=TE2-TF1 

TEEROR=ABS( TEDIFF)/T?:1 

IF (TEBROR-C ,n l)7n,3f ,,3'. 

30 IF (XCONVl-XCONV2)5t,, 60,40 
4i K=K+1 

TELOW=TELOW 

TEHIGH=TE2 

TeRANG=TEHI6H-TEU0W 

TELOCT=Ci, 382C«*TERANG 

TE2=Tei 

TI2=TU 

XI2=XI1 

XC0NV2=XC0NV1 

Tei=TFLOW+,T£LOCT 

CALL MATENG(XE,Tai,P,YF,ACP,BCP,CCP,F,XP,OP,XIl,TU,VI,VEt 
1TIM0AY,H2SPPM) 

XC0NV1=XE-XII 
GO TO 20 
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5t,.; K=K + 1 

TEL0W=TE1 

TSHIGH=TEHIGH 

T5RANG=TEHIGH-TEtjOW 

TELOCT=n» 3e2t »TERANG 

TE1=TE2 

TI1=TI2 

XI1=XI2 

XCQNVl = XC0f^lV2 

T52=TEHI6H»TEL0CT 

CALL MATEN3(XH,TB2, P,YF,ACP,8CP,CCP,FtXP,DP,XI2 ,TI2»VI» VEt 
1TIMDAY,H2SPPM) 

XCQNV?.= Xe-XI 2 
GO TO 2C 
6( K=K+1 

TFLnW=T£l 
TEHIGH=TE2 
GO TO 16 

TC I F ( XCON VI -XCONV 2 ) 90 1 90 ,.8'? 

8'; XCOPTM^XCONVl 
TeaPTM=T6I. 

TI0PTM=TI1 
XI0PTM=XI1 
GO TO im 
90: XC0PTM=XC0NV2 
TE0PTM=TE2 
TiaPTM=TI2 
XI0PTM=XI2 
GO TO 1CC5 
It i: CONTINUE 
16L RETURN 
END 

subroutine RATMCH{XIN,TIN,P, YF»XP,DP,RATIN,TENltTIM0AY,H2SPPM) 

— RATE-MATCHING PROGRAM SATISFIES FOLLOWING CONDITIONS BETWEEN TWO 
— - BEDS (N-1 )AND N. 

— - CONDN.l RATE AT EXIT OF BED { N-l )=RATE AT INLET TO BED N. 

— - C0NDN.2 XE AT EXIT OF BED (N“1)=XI AT INLET TO BED N. 

— C0NDN.3 TEMP.AT EXIT OF BED (N-l) IS GREATER THAN TEMP. AT INLET TO BED 
DIMENSION YF{ 10) 

XMEAN=XIN 

CALL RATEQN(XIN,lIN,Pf YF,XP,DP,XeQLIN,RATIN,XMEAN,TlMDAY H2SPPM) 

— N1=NL = BED.(N~I) 

■— N2=BED N 
XENL=XIN 
DELTAT=2?5.P 
I'O K=C 

TENL=TIM+DELTAT 
15 K=K+1 

CALL RATEQN {XENL,'TENLtP,YF,XP,DP,XEQLNL»RATENL, 

:XMEAN,TIM0AY,H2SRPM) 
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2:; RATDIF = RATENL-RAT!IN 

RTeRQR=A8S( RATDrF)/RATIM 
I F ( ftTER.OR.-!>e<^«^l ) 50, 25, 25 
25 IF (RATDI F) 3.N50,'35 
3' DELTAT = DELTAT«! o25r 
IF(K-l)ir»,l(.,4J. 

35 TEOLD=TENL 
4* TENL=TEOLD+DELTAT 
GO TO 15 

C*** FOLLOWIMG COMDNS« CORRESPOND TO EGUAL RATES AT EXIT OF BtO N1 

C«*» AND '=,NTRANCF OF BFD N2e 
5i:; RATIN2=RATIN 
RATEN!=RATENL 
XEN1=XENL 
TENl=TeNL 
RETURN 
END 

FUNCTION TEQ(YF,X) 

C- — THIS FN, FINDS TEQ (cQUI LoTEMPo ) IN DEGREE K AT A GIVEN DEGREE OF 

C CONVERSION X AND FEED COMPOSITION YFCMOLE FRACN.) 

DIMENSION YF(V ) 

XEQUIL=X 

YCONV=YF( 2)*XEQUIL 

p.. EQUILK = i { YF (3)+YC0NV)*{ YF (1) +YCONV) }/( { YF (2 )”YCONV ) *(YF {7)»YC0NV ) ) 

f c FINDING TEMPoCORRESPONOING TO EQUI LK( ACTUAL VALUE OF EQUI L, CONST* ) 

: C BY TRIAL AMD pRRCR SEARCH* 

TMAX=9uo.r'+273o 
DELTAT = 2D‘'.e(: 
lo KI=M 

TEGGrS=TM AX-DtiLTAT 
IS KI=KI+l 

EQKGES=EQK{ TEQGES ) 

C EQKGES IS VALUr OF EGUI L« CONST* AT A GUESSED TEMP* TEQGES* 

2C 6QKDIF=EQKGES-EGUILK 

EKEROR = ABS{ eCKD IF )/EQUILK 
IF (EKEROR*' *'0.1)5’ ,25,25 
25 IF {EQKDIF)35,5<5,'3n 
3€ DeLTAT=DELTAT*C«25U 
IF (KI-1) 10,10,4C 
35 TeQOLD=TEQGES 
4C TEQGES=TEQ3LD-D£UTAT 
GO TO 15 
5f: TEQ=TEQGES 
RETURN 
END 

SUBROUTINE MATEN.G ( X I , TI ,P , YF , A , B ,C ,F , XP ,DP ,XE,TE, VE , VI , 
ITIM0AY,H2SPPM ) 

C*** SIMULATION OF SHIFT REACTORI FINDING CONV.AND TEMP* FROM GIVEN VOLUME) 
PERFORMS MATERIAU-ENSRGY BALANCE*V IS INDEPENDENT VARIABLE* 

THIS SUBROUTINE SOLVES TWO SIMULTANEOUS FIRST-ORDER DIFFERENTIAL 



o n o 


95 


EQNS, USIM3 RUNGE-KUTTA METHOD. FINDS X(DEGREE OF CONV. OF CO) 

AND T(TEMPo) AS A FN.DF ViCATALYST VOL.) BETWEEN VI=fi AND VOLUME { N8 ) 
INITIAL CONDITIONS VI,XI,TI ARE GIVEN. 

DIMENSION A( ■7),B( ■^),YF{lL) 

dimension x{it),it(i f.' ) ,v ( 1*”-: ) 

DIMFMSION K’rCLdO- ) tSAT'^ {?..!■'> ) ,TTRY 
K=;- 

MINTVL=5 
G5 K=K+1 

NINTVL=NINTVL + E 

1’ OFLTAV=(VF"VI )/F.LOAT(NINTVL) 

JNC=0+l 

JNINT=NINTVL+I 

X(JNO)=XI 

TUNr])=TI 

V(JNO)=VI 

DO 20 JM= JNfl, JNINT 
XMEAM=X( JN) 

N=JN"I. 

XJNi=X{JN) 

TJN1=T(JN) 

CALL RATEQN(XJNlvTJNl,PtYF,XP,DP,X£QJNlfRATJNl, 

3XMEAN,TIMDAY,H?.SRPM) 

ALl=Or:LTAV»DTBYC)?(XJNl,TJNl,P,YF,A,B,C)*DXBYDVlRATJNl,F,YF) 
AK1=DFLTAV*DXBYDV(R ATJNl, F,YF) 

XeQL{4N)=XECJMl 
RATE( JN)=RATJNl 
XJN2=X(JN) + AKl/2. : 

TJN2=T( JN ) + ALl/2.C 

CALL RATEQN(XJN2,TJN2,P,YF,XP,0P,XEQJN2,RATJN2f 
1XMFAN,TIMDAY,H2SRPM) 

AL2=DELTAV*DT8YDX(XJN2»TJN2, P,YF,A,8,C)*DXBYDVtRATJN2,FtYFy 
AK2 = DPLTAV»0X8YDV(RATJN2tF,.YF) 

XJN3 = X( JN)+AK2/2.''. 

TJN3=T( JN )+AL2/2.r 

CALL RATEQN(XJN3,iTJN3,P,YF,XP,DP,XtQJN3,RATJN3, 

IXMEANtTIMDAY, H2SRPM ) 

AL3=DELTAV»DTBYDX(XJN3,TJN3t P»YFfAf B,C)*DXBY0V{RATJN3,FtYF) 
AK3=DELTAV*DXBYDV(RATJN3,F,.YF) 

XJN4=X( JN ) + AK3 
TJN4=T( JN)4-AL3 

CALL RATeQN(XJN4,'TJN4,P,YF,XP»DP,XeQJN4,RATJM4, 

IXMEAN,TIMDAY,H2SRPM) 

AL4=DELTAV*DT8YDX(X JM4tTJN4, P,YF ,A tB , C ) *DXBYDV { RAT JN4t F , YF ) 
AK4=DELTAV»DXBYDV(RATJN4»F,.YF) 

DELTAX=(.AKI+2«u»AK 2+2.C*AK3+AK4) /6.fi 
DELTAT={AL1+2.0«/!L2 + 2.':.'»AL3 + AL4) /6.C 
T{ JN+1)=T(JN)+DEUTAT 
V( JN + 1 )=>V{JN)+DE.UTAV 
X( dN + 1 )=X(JM)+0 -.UTAX 



2«^ CONTINUE 

TTRY(K)=T (JNINT) 

IF (K- 1)35, 3 5, At; 

AH EftOR=ABS{ TTRY{K HTTRY{ K-1 ) ) 

IF( (BROR/TTRY(K) I-« ; 1) '3* 1 35 , 35 

5c TF=TTRY{.K) 

XH=X{JNIMTy 
XECLE=XFQL( JNINT 1 
RATEE=RATE( JNINTI 
XEQLI=XFQL( JMO) 

RAT-I=RATP{ JMQ) 

RKPRMT = 1«»;. 

IFtRKPRNT )9( , 55,9 
55 PRINT 6f 

6. FORMAT{///5X,«RUNGF“KUTTA INTEGRATION*) 

PRINT 7 5, XI , X FQLI ,T I , V I ,R ATgl 

73 FORMAK 5X,»XI=*,P7,4,5X,*Xr;QLI=* ,F7. A , 5X, *TI =* , F7» 2 ,5X, 
l*VI=s*,.F7» 3, 5X,*RATEI=*,E15«3/) 

PRINT 8t',XE,X=QLE,TE,VB,RATEE 
8G FORMAT{5X,»XF=*,F7,4,5X,*Xr:QLE=*,F7.A,5X^*TE=*, 
iF7,2t 5X,*VE=*,Fl.*3,.5X,*RATEE = «,E15o3///) 

9r CONTINUE 
RETURN 
END 

FUNCTION DT3YDX(X,T,P,YF,A,BrC) 

DIMENSION YF(il;),*A( 7),B( 7)tC{ 7 ) , Y( 100) ,CPO ( 1€3 
DATA U,V, W, AKl, AK2, TO, AD, 8D,CD/0,63E14,.-5.56,1» 29,3«g33Et4 
l-0,93 539E-f'3,6C'!.'n,-^''..2 .1,5. A93E-C?3,-2«7177E-06/ 

CALL MASS (I ,YF, K,'Y) 

cposuM=r.4t 

DO If! 1=1,7 

CPO( I) = A{ iy+B{I)<(T+C(I )*(T**2.C ) 

CPOSUM=CPOSUM+Y (I )*CPO{ I) 

If CONTINUE 

CPRSUM=CP0SUM+Y(7 )*U*(P«*W)*{T**V) 

HRTOP=-AK 1» EXP ( AK2*P ) 

Cl = ( { AD*TO) +{ BD«LT0**2.t ) ) /2*f, + ( CD»{T0**3.©) ) /3»t ) 
C2=U/(V+1.0) 

C3=(U/( V+l.C) )*{10*«(V+'1.€)) 

HRTP = HRTOP+;( ( AD«1 ) + i BD* ( T**2*0) / 2.f' ) + { C0» I T»*3. 0/3.0) 
1-{C2*( P**W) *( T**LV+1.P) )) )+'C3*( P«*V»)--C1 
DTBYDX=t-HRTP*YFi2) )/CPRSUM 
RETURN 
END 

FUNCTION DX8YDV(RATE,F,YF) 

DX8YDV=DX/DV 
DIMENSION YF(ir) 

DX6Y0V=RATt /I F*YR (2 ) ) 

RETURN 

END 
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SUeaOUTINF aATFQ^(XtT,P,YF,XP,DP,XSQL,RATe»X^*EAN,TIM 0 AY,H 2 SPPM) 
RATEQN CALCULATES PSEUDO FIRST-ORDER REACTION RATE 0BTAIM5D 
C*** BY SIMPLIFYING SECOND ORDER RATE E XPRESSI ON ( L AUPICHLER» MOE) » 
e«** BASED ON RUTHVENS METHOD FOR FCI HT SHIFT CATALYST 
C«»* RATE EQN IS DFVEliOPFD FOR PRESSURE RANGE 6-26 KG/SQ.CM.ABS 
C«** UNITS — TIMDAY=TiMF IN DAYS! ON STREAM) — (GREATER THAN I DAY) 

C«** UNITS— H2SPPM = CCNCNoOF H2S IN PFMoCDRY BASIS) 

C«*» UNITS— TFMPFRATLRF=t DEGPEF K, 

C^|** UNITS— PRFSSURE=?P KG/SQoCM. ABS 

C*** DP = cFFeCTIV- PARTICL’* DIa'IN CM« =6*VP/AP 

C**» RK STANDS FOR RATE CONST*, K 

C * * «»«**«• IS *•»■»**** ««*#*#**#*#»*#*«.« 

DI Ms=NS ION Y F { If ) ,<YMEAN( ini ) 

C — - ZCORCN = CORR,ECTIGN FACTOR FOR RATE CONSTANT. 

ZCORCN=l.f 

R=1.9B7 

Z0=4.22I"''B 

ZTIMF=TIMDAY*#(“0«C79) 

IF (H2SPPM-2«7)15,16,16 
15 ZSULFR = l.r 

i GO TO 2e 

' 1C ZSULFR=((-D.276)<lALOG10(H2SPPM) ) +1.127 
20 CONTINUE 

^ FA=13.43r» 

■ ZPRESS=t5XP (-FA/P)) /P 

cA = (9.225rC3)»P 
ZTeMP=FxP(-FA/( R+tT) ) 

RKV=Z0*ZTIME*ZPRBSS*ZTEMP*ZSULFR«ZC0RCN 

CIFFPl='\r;43 

0IFF=(DIFFPl/P)*t(T/673«0)»*l,5) 

THMOD = ( DP/2.T')*SQRT(PKV/DIFF) 

IF (THMOD)35t25r35 
25 £FF*l.(i 
GO TO 40 
35 CONTINUE 

EFF»( 3.ft/THMOO)*t(l.r/TANH(THMOO))-{l.t?/THMOO) ) 

4C! CONTINUE 
V0ID=«).4 

RK0=RKV*EFF*(1.C-V0I0) 

RKAi=273.0*3 6PC««';-H(P/T)*RK0 
YMEAN(7)=.YF(7)-YP(2)*XMEAN 
RATE=(RKA*YF( 2) )^^ (XEQ{ YFtT)-X)*Y#'EAN(7) 

XEQL=XEQl YF,T) 

RETURN 

END 

FUNCTION XEC(YF,1) . 

C THIS FN, CALCULATES XEQCDEGR6E OF CONVERSION OF CO) AT A GIVEN 

C TEMP. T DEGREE K AND FEED COMPOSITION YFtMOLS FRACN.) 

DIMLHSION YF( 10 ) 

AI=EQK(T)-l.r 


* # 



98 


B1={-£QK( T)«YF(2) )“(^QK{T)»YF(7) )“YF(3)-YF(1.) 

Cl={EQKtT )»YF{2) ^^YF {?))”{ YFt 3 )*YF(1) ) 

YC0NV2 = (-Bl-'SQRTt{Bl#81)»{4e; •»A1«C1) ) )/(2.y*Al) 

YCQFQ2=YF(2 )“YCCNV2 
IF(YCaSQ2 )5i', 25,25 
25 YC0£Q=YC0EQ2 
GO TO 3r 

S* YCONVl={-Bl+SQRT(.{Bl<fBl)«»{4.r#Al«Ci)n/(2,r,*Al) 

YC0L'C1=YF (2 )-YCCNVl 
IF(YCnEQl)7^ ,2*; ,2.' 

2t YCOFQ=YCOeQl 

3! XrG = (YF{2 )»YCaFG)/YF(2) 

XEGL=XEQ 
FQUILK=r_QK( T) 

GO TO 3f 

7: PRIM 75, YC0'^G1,YC05G2 

75 FOftMATC 5X,*YC0EG1=«,F7»4, 5X,*YCOEQ2=*,F7,4,5X, 

1*TFIS IS NOT POSSIBLE*) 

8C RETURN 
cNC 

FUNCTION EQK(T) 

C THIS FNo CALCULATES EQUILIBRIUM CONSTANT EQK FOR WATER GAS SHIFT 

C REACTION AT A GIVEN TEMPERATURE T CgGREE K. 

DATA R, A0,B0, CD, HOO, CON SI , SDO,CONS2/l, 987,-0* 201 , 5» 493E-D3 , 
l~2«7l77E-r*6,-9B38o0,160*22C,-10»C42,.0,372/ 

HDT=HDO + ( (AD*T) +{.BD»T*T/2*0) +(CD«{T**3«*l)/3*e) )-CQNSl 
S0T=S0Q4( {AD*ALCG{T ) )+(BD*T)+(CD*T*T/2.0) )”C0NS2 
FDT=HOT-t T*SDT) 

SQK=EXP{-FDT/(R,*7) ) 

RETURN 

END 

SUBROUTINE MASS { NCI , YF ,.XT, Y) 

C*»* STOICHIOMETRIC RBLATIONS 
DIMENSION YFdJ ),'Y{l<r:) 

YC0NVT^YF(2 )*XT 
Y(NC1)*YF{1 ) + YCONVT 
Y(NCl + l)=<YF(2)-YaONVT 
Y(NC1 + 2)=!YF(3)+YOONVT 
Y( NCl+3)=YF{4) 

Y(NC1+4)=YF(5) 

Y(NC1+5)=*YF(6) 

Y( NC1+6)=YF(7)“YQ0NVT 
RETURN 

END ' ■ . . 


29 »0C 

60 

3 

53»2r.- 

46 '*6^’ 

3* 20 

500 e f-f 



4* 35'^ 


8 ,7.^0 

45* 

27'','’ 


i C’- 

1 




99 


BErBKSlGES 

1. Denbigli, K. G. and Turner, J.C.R. ; '’Chemical Reactor 
Theory” Cambridge University Press, (1971), 

2. Horn, P. and Euchler, L. ; Chem.Ing. Techn. 1(1959). 

3. Bilous, 0. and iimundson, N. R. ; Chem. Sngng. Sci., 5, 

81, 115 (1956). 

4. Kramer, H. and Westerterp, K.R. j "Elements of Chemical 
Reactor Design and Operation” Netherlands University 
Press, Amsterdam (1963), 

5. Laupichler, E. G. ; Ind. Eng. Chem., 30, 578 (1938). 

6. Padovani, C. and lotteri, A. j J, Soc, Chem. Ind., 

391 T(1937). 

7. Kulkova, N.V. and Temkin, M.I,, Zhur. Eiz. Khim. 

695 (1949); Chem. Abstrs. , 43, 7308h (1949). 

8. Atwood, K. , Arnold, M.R. and Appel, E, G. , Ind. Eng. Chem, , 
42, 1600(1950). 

9. Mars, Cliea. Engng, Sci,, H, 375 (1961). 

10. Popov, B.I., Zhur. Eiz. Khim. , 31., 1033 (I95l). 

11. Barkley, I.W. , Corrigan, T,E. , Wainwright, H.¥, and 
Sands, A.E. , Ind. Eng. Chem., 1066 (1952). 

12. Kodama, S, , Eukui, K. , Tame, T. ; and Kinoshita, M. , 
Catalysts (Japan), 8, 50 (1952). 

13. Kodama, S, , Mazume , A. , Eukuba, K. and Eukui, K. , Bull. 
Chem. Soc, Japan, 318 (1955). 

14. Bortolini, P. , Chem. Engng. Sci., 135 (1958). 

15. Stelling, 0., and Krusenstierna, 0. Von, Acta Chem. Sc and, , 
12, 1095 (1958). 



100 


Slrrlllov, I.P., ZM..-Iektool. Inst.. 

Ho.S, 45 ( 1956 ); Chem.^Abstr. 53, I1969 (1959). 

17. ■^'ti'oshciienlco. Y T 

• • » ' -- ibr. 3 . , 21 iur. Priklad, Khim. 

12. 997 ( 1959 ); Ohem. Abstr. 15735 (1959)_ 

18 . Hulburt, H.M., ana Tasan, C.D.S., A.I.Ob.B. Journal, 7 
143 ( 1961 ). 

19 . Boblbro, H., Acta Ohem. Seand. , _15, 502 (1961). 

20 . Atroabchenlco, y.i., Zaaorin, ' a.p. 

Alnetijoa 1 Katallz, 5, 605 ( 1962). English Translation 
in Kinetics and Catalysis, 3 529 (1962). 

21 . Ooodridge, P. , and Quasi, H.A., Trans. Inst. Ohem. Sngrs. , 
45 , T 274 ( 1967 ). 

22 . Moe, J.M., Ohem.Eng. Progr. , lfo. 3,33 (1962). 

23 . Hougen, O.A., Watson, K.M., and Ragats, R.A., '.Ohemioal 
Ptocess Principles.., Part I, John Wiley Sons, Inc. (1954). 

24 . Hilsenrath, J. , Beckett, O.W., Benedict, w. S. , i^no, I., 
floge, H.J., Basi, J.p., Uuttall, R.I., loulofcian, T. S. , 

Wolley, H.w. , "Tables of Thermodynamic & Transport 
Properties" Pergamon press, Oxford (i960). 

25 . Huthyen, D.M. , Ohem. Eng. Scl. 759 (i9g8)_ 

26 . Ruthven, D.M., Oand. J. Ohem. Eng., £7, 327 (1969). 

27. Thiele, E.W., Ind. Eng. Ohem., £1^, gig {1939) 

28 . Hoogschagen, J. , md. Eng. Ohem., 47, 906 (1955). 

29. Satterfield, O.H. and Sherwood, T.E., ..The Role of 
Diffusion In Catalysis.', Addison-Wesley (1963). 



101 


30. Chandra, m. , singh S S 

, •’ B., Sinha, A.z and 

Sinha, If, 1C 

» ape r presented at ’’The a 

MeetiiiP- Tr,^- Annual General 

S Indian Institute of C^’enjoai p • 
at TTrp Tr -Smical Engineers" held 

at III-Kanpur (Peh. 1972). 

31 . Mahapatra, h. , Bay, jf se„ ^ 

y, n., ben, B. and Sen, s.p 'Pa^hr.^T 
8, 211 (197^)^ -iechnology, 

ling, A.p, and Wan, S. W Cham i? 

33 . Banerjee 3 oha rt ' (1969). 

’ °1iandra, M. and Shosal, 3 R „ 
at "The 2ith A J-aper presente 

-^ae 24th Annual General Meetiu.^ or r 
of GhPm- ^ i^eeting of Indian Institute 

iiemical Engineers" held at ttti / 

34. LavroT, u t Ror v. Hl-Sanpur (Jeb. 1972 )., 

•’ aiiPPOva. y.i.. 

Ihermo dynamics of Ga si float t 

Per.. ^^safication and Gas-Synthesis Reactions" 

Pergamon Press (1963). -^^i^ions 

35. Smith, J.M., "Chemioal Bngiaeeria, r- r 

-^^(giueering Kinetics". McGr.t.t w,-t 7 
Book Company (1956), ~ ^ 

'^Hemical Bn^ineera.. 

iicbraw-Hili, New York ( 1962 ). 

37 . Aris, R. , "Optimal Design of Chemi i d 

Press, New York (1961). factors", Academic 

38 . Horn, 7 . P. , Z. Eiektrochem. , . 295 (l 96 l). 

a.-al, M. D, , ard Dranoff, j s ind F n 
^ 3 ign Develop., 5, 129 (1966). ^ocess 

W. Beveridge, 0,3. and Schechter Rs un t- • 

“Py ana Practice" Moeraw-Hill, Rew York (1970). 



